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ABSTRACT 

Search  of  an  unknown  space  by  a  physical  agent  (such  as  an  autonomous  vehicle)  is 
unique  in  search  as  the  customarily  most  important  goal  (to  reduce  the  computation  time 
required  to  obtain  the  shortest  distance)  is  not  as  important  as  minimal  movement.  There  is 
a  real-time  aspect  since  the  agent  is  actually  moving;  using  energy  each  step  of  the  way. 
Having  limited  energy  resources  and  knowledge  of  the  terrain  (only  adjacent  nodes),  the 
key  factor  for  the  physical  agent's  search  algorithm  is  reduction  of  steps.  Hence,  any 
heuristic  that  can  help  keep  step  count  to  a  minimum  must  be  considered.  Korf  and  Shing 
addressed  this  issue  in  separate  works.  Both  made  use  of  known  information  about  the 
frontier  node's  distance  from  the  current  node  in  addition  to  a  heuristic  estimating  the 
distance  from  goal. 

In  this  thesis,  we  present  a  simple  genetics-based  method  to  produce  adaptive,  efficient 
multi-heuristic  search  strategies  for  the  real-time  problem.  Extensive  empirical  study 
shows  that  this  approach  produced  search  strategies  with  much  better  performance  over 
existing  search  algorithms  for  most  terrain  types.  The  methodologies  used  to  develope 
these  improved  strategies  for  our  specific  case,  are  also  applicable  to  a  multitude  of  real- 
time search/optimization  problems  in  the  general  case. 


in 


CI 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 1 

II.  PROBLEM  MODEL 3 

A.  TERRAIN  MATRIX 3 

B.  DENSITY  MATRIX 3 

C.  TERRAINS  USED 4 

1.  Central  mountain 4 

2.  Single  Left  Ridge 5 

3.  Single  Right  Ridge 5 

4.  Double  Ridge 6 

5.  Single  Left  Plateau 7 

6.  Single  Left  Plateau  with  Ridges 7 

7.  Random  Terrains 8 

ffl.     BACKGROUND  WITH  DISCUSSION  PERTINENT  TO  MODEL 9 

A.  KNOWN  TERRAIN  SEARCH 9 

1.       A*  Search 9 

B.  UNKNOWN  TERRAIN  SEARCH  BY  A  PHYSICAL  AGENT 9 

1.  Hill-climb  Search 10 

2.  Real-time-A*  Search 10 

3.  Best-first  Search 10 

4.  Persistence  Search 11 

C.  GENETIC  ALGORITHMS 11 

IV.  FACTORS  RELEVANT  TO  SEARCH 14 

A.  DISTANCE  FROM  START 14 

B.  DISTANCE  FROM  CURRENT 14 

C.  DISTANCE  FROM  GOAL 14 

D.  CROWDING 14 

1.  Crowding  sides 15 

2.  Crowding  diagonals 15 

E.  MOVE  AWAY  FACTOR 15 

F.  MOMENTUM 15 

V.  PROGRAM  DEVELOPMENT 16 

A.  DATA  STRUCTURES 16 

1.  Node  structure 16 

2.  Population  structure 16 

3.  Frontier  Heap 17 

4.  Mate  Heap 17 

B.  SEARCH  ALGORITHMS 17 

1.  A*  Search 17 

2.  Hill-Climb  Search 17 

3.  Real-Time-A*  Search 20 

4.  Best-First  Search 20 


IV 


"  '^^\  i-iorvA\r\  r 


NAVAL  POSTGRADUATE  SCHr.ni 
MONTEREY  CA    939«  5101 


5.  Persistence  Search 20 

6.  Multi_heuristic  Search 24 

C.      GENETIC  ALGORITHM 26 

VI.  TRAINING 28 

VII.  TESTING 30 

VIH.  EXPERIMENTAL  RESULTS 31 

A.  NATURAL  TERRAINS 31 

1.  Central  Mountain 31 

2.  Single  Left  Ridge 31 

3.  Single  Right  Ridge 33 

4.  Double  Ridge 34 

5.  Single  Left  Plateau 35 

6.  Single  Left  Plateau  With  Ridges 36 

7.  General  Comment 36 

B.  RANDOM  TERRAINS 37 

1.  Random  One  /  Random  Two  /  Random  Three 37 

2.  Random  Four 39 

3.  General  Comments 40 

C.  GENERAL  COMMENTS 40 

IX.     CONCLUSIONS 41 

LIST  OF  REFERENCES 42 

APPENDIX  A 

TERRAIN  DEVELOPMENT  FROM  A  SAMPLE  DENSITY  MATRIX 43 

APPENDIX  B 

RANDOM  TERRAIN  DENSITY  MATICES 45 

APPENDIX  C 

SEARCH  STRATEGY  COMPUTATIONAL  TIME 47 

APPENDIX  D 

PROGRAM  C  CODE 48 

INITIAL  DISTRIBUTION  LIST  124 


LIST  OF  FIGURES 

Figure  1  Central  Mountain  Terrain  4 

Figure  2  Single  Left  Ridge  Terrain  5 

Figure  3  Single  Right  Ridge 6 

Figure  4  Double  Ridge  Terrain  6 

Figure  5  Single  Left  Plateau  Terrain 7 

Figure  6  Single  Left  Plateau  with  Ridges  Terrain 8 

Figure  7  Individual  Chromosome  Structure 16 

Figure  8  A*  Search  Algorithm  18 

Figure  9  Hill  Climb  Search  Algorithm 19 

Figure  10  Real  Time  A*  Search  Algorithm  21 

Figure  1 1  Best  First  Search  22 

Figure  12  Persistence  Search  Algorithm 23 

Figure  13  Multi  Heuristic  Search  Algorithm  25 

Figure  14  Allele  Crossover  Example  26 

Figure  15  Bit  Crossover  Example  26 

Figure  16  Genetic  Algorithm  27 

Figure  17  Training  Algorithm 29 

Figure  18  Central  Mountain  Results  32 

Figure  19  Single  Left  Ridge  Results  32 

Figure  20  Single  Right  Ridge  Results 33 

Figure  21  Double  Right  Left  Ridge  Results 34 

Figure  22  Single  Left  Plateau  Results  35 

Figure  23  Single  Left  Plateau  With  Ridges  Results 36 

Figure  24  Random  One  Results  37 

Figure  25  Random  Two  Results 38 

Figure  26  Random  Three  Results 38 

Figure  27  Random  Four  Results 39 

Figure  28  Sample  Density  Matrix  43 

Figure  29  One  of  Many  Possible  Resultant  Terrains  44 

Figure  30  Random  Terrain  One  45 

Figure  31  Random  Terrain  Two 45 

Figure  32  Random  Terrain  Three 46 

Figure  33  Random  Terrain  Four 46 


VI 


I.  INTRODUCTION 

Search  of  an  unknown  space  by  a  physical  agent  (such  as  an  autonomous  vehicle)  is 
unique  in  search  as  the  customarily  most  important  goal  (to  reduce  the  computation  time 
required  to  obtain  the  shortest  distance)  is  not  as  important  as  minimal  movement.  There  is 
a  real-time  aspect  since  the  agent  is  actually  moving;  having  limited  time  to  determine  its 
next  move  and  using  energy  each  step  of  the  way.  This  is  in  contrast  to  the  traditional 
problem  of  search  of  known  space  for  the  shortest  path  which  can  be  efficiently 
accomplished  by  A*  search  with  a  good  heuristic  estimating  the  distance  to  goal.  The  path 
is  found  without  any  movement.  Although  factors  other  than  the  actual  distance  from  start 
and  estimated  distance  from  goal  could  reduce  the  number  of  nodes  examined  in  the 
traditional  problem,  these  factors  usually  increase  the  computational  cost  per  node 
examined  and  produce  paths  that  are  longer  than  the  shortest  path  which  makes  additional 
heuristics  undesirable.  Such  is  not  the  case  for  the  real-time  problem. 

The  physical  agent  traversing  a  terrain  in  the  real-time  problem  knows  only  its  current 
position,  the  goal's  position,  and  whether  adjacent  and  previously  adjacent  nodes  are 
passable  or  not.  It  learns  about  the  terrain  only  as  it  moves  from  node  to  node  examining 
all  nodes  adjacent.  Information  about  past  nodes,  visited  or  adjacent,  can  be  stored  to  build 
up  its  knowledge  base.  Computational  time  to  determine  the  next  move  is  important,  as 
stopping  to  compute  before  each  move  is  undesirable.  On  the  other  hand,  insufficient 
computations  can  result  in  unnecessary  steps  and  wasted  energy. 

Having  limited  energy  resources  and  knowledge  of  the  terrain  the  key  factor  in  the 
physical  agent's  search  is  the  reduction  of  physical  steps.  In  [Pa89],  Papadimitrion  and 
Yannakakis  showed  that  the  computational  problem  of  deriving  optimal  search  strategies 
for  the  real-time  problem  is  PSPACE -complete.  Hence,  any  heuristic  that  can  help  keep 
step  count  to  a  minimum  must  be  considered.  Korf  [Ko90]  studied  this  problem  and 
developed  the  real-time-A*  search,  which  uses  the  physical  agent's  distance  from  the  node 
(g  (n) )  in  addition  to  the  distance  from  goal  heuristic  (h  (n) )  to  determine  the  best  next 


move  by  minimizing  the  objective  function  f(n)  =  g  (n)  +  h  (n)  for  every  adjacent  node 
n.  Shing  and  Mayer  [Sh91]  developed  persistence  search  which  included  a  persistence 
factor  (pf  =  0  to  1)  to  bias  the  distance  from  current.  The  next  move  is  determined  by 
minimizing  the  objective  function  f(n)  =  pfx  g  (n)  +  h  (n)  for  every  frontier  node  n. 
Experimental  results  led  to  the  conclusion  that  the  pf  could  be  adjusted  to  optimize  search 
depending  on  terrain  type  and  the  density  of  obstacles.  Details  of  these  search  strategies  are 
in  Chapter  in. 

Extending  on  these  works,  we  believe  a  combination  of  additional  heuristics  can  be 
beneficial  in  minimizing  physical  agent  steps.  As  the  number  of  heuristics  increases,  it  is 
essential  to  have  some  efficient  means  of  assigning  bias  adjustments  to  various  heuristics 
to  optimize  f( n)  for  different  terrain  types  and  densities  of  obstacles.  If  the  combinatorial 
explosion  required  to  produce  all  possible  combinations  of  heuristics  is  not  intractable,  the 
required  testing  of  each  to  select  a  best  makes  this  means  computationally  prohibitive. 
Since  enumeration  is  probably  not  possible,  some  random  means  of  attaining  the  best 
combination  seems  to  be  the  most  plausible.  DeJong  [De75]  made  clear  the  advantages  of 
genetic  algorithms  over  purely  random  selection. 

In  this  thesis,  we  present  a  simple  genetic  algorithm  based  method  to  produce  adaptive, 
efficient  multi-heuristic  search  strategies  for  the  real-time  problem.  Extensive  empirical 
study  showed  that  this  approach  produced  search  strategies  with  much  better  performance 
(reduced  number  of  steps  without  prohibitive  computation  time)  over  existing  search 
algorithms  for  most  terrain  types.  The  methodologies  used  to  develope  these  improved 
strategies  for  our  specific  case,  are  also  applicable  to  a  multitude  of  real-time  search/ 
optimization  problems  in  the  general  case. 


H.  PROBLEM  MODEL 

A.  TERRAIN  MATRIX 

To  best  demonstrate  the  effectiveness  of  the  multi-heuristic  search  strategies  produced 
by  a  genetic  algorithm,  we  chose  to  apply  the  strategies  to  random  obstacle  distributions  in 
the  form  of  a  two-dimensional  64x64  grid  of  squares  (nodes).  Nodes  can  be  either  free  or 
obstacles,  movement  can  be  in  eight  directions  through  free  spaces  only.  A  perimeter 
surrounding  this  grid  is  a  solid  row/column  of  obstacles.  The  distance  from  a  node  to  its 

horizontal/vertical  neighbor  is  1.0;  to  its  diagonal  neighbor  is  Jl.  The  total  distance 
traveled  from  start  to  goal  according  to  any  search  scheme  is  the  sum  of  each  of  these 
individual  steps.  The  effectiveness  (fitness)  of  a  specific  search  scheme  is  the  ratio  of  the 
shortest  path  length  from  start  to  goal  divided  by  the  distance  traveled.  Given  as  a 
percentage,  100  is  the  best  possible;  meaning  the  distance  traveled  is  equivalent  to  the 
shortest  path.  Specific  nodes  of  the  grid  are  be  identified  by  Cartesian  coordinates  with  the 
left  border  column  being  the  v  axis  and  the  bottom  border  row  being  the  x  axis.  The  lowest 
left  node  is  (1,1);  the  top  right  is  (64,64). 

B.  DENSITY  MATRIX 

The  64  by  64  search  space  grid  is  divided  into  16x16  density  blocks,  each  containing 
4x4  nodes  and  having  a  specified  block  density.  Block  densities  range  from  0-15.  A  block 
density  of  9  means  that,  on  average,  nine  of  the  block's  16  nodes  will  be  an  obstacle  (chosen 
at  random).  These  density  blocks  are  numbered  from  (0,0)  to  (15,15)  where  (0,0)  is  the 
bottom  left  and  (15,15)  is  the  top  right.  Start  and  Goal  positions  are  specified  by  density 
blocks.  Most  of  the  block  density  distributions  used  will  have  a  start  block  of  (2,2)  and  a 
goal  block  of  (13,13).  The  specific  start/goal  node  is  located  randomly  in  that  block.  See 
Appendix  A. 


C.     TERRAINS  USED 

There  are  ten  different  density  distributions  that  were  used  for  training  and  testing.  The 
block  densities,  once  set,  remain  unchanged  from  the  start  of  training  through  testing. 
Although  the  block  densities  remain  constant,  actual  obstacle  placement  is  determined 
stochastically  and  changes  from  run  to  run.  The  point  is  to  investigate  the  adaptability  of 
genetic  algorithm  to  produce  the  best  strategy  to  direct  the  search  through  terrains  where 
the  general  density  distribution  is  known  but  actual  obstacle  placement  is  not.  The  first  six 
terrains  are  considered  natural  terrains  since  they  closely  resemble  actual  topological 
conditions.  The  start  density  block  is  always  (2,2)  unless  otherwise  stated.  The  goal  density 
block  is  always  (13,13)  unless  otherwise  stated. 

1.     Central  mountain 

The  highest  density,  15  (denoted  as  f  in  Figure  1),  is  in  the  center  with  a  gradual 
decrease  towards  the  lowest  density,  1,  on  the  outer  edge.  Figure  1  shows  the  density 
distribution  of  the  terrain  in  hexadecimal.  Transit  from  start  to  goal  requires  a  search 
scheme  to  find  the  most  efficient  way  around  the  mountain. 
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Figure  1  :    Central  Mountain  Terrain 


2.     Single  Left  Ridge 

This  terrain  has  a  high  density  (15)  ridge  starting  from  left  center  moving 
horizontally  out  past  the  grid's  midpoint.  There  is  a  gradual  decrease  in  density  down  to  2 
as  the  distance  increases  from  ridge  center  (Figure  2).  Transit  through  the  ridge  is  not 
possible. 
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Figure  2  :  Single  Left  Ridge  Terrain 


3.     Single  Right  Ridge 

This  terrain  is  similar  to  the  Single  Left  Ridge  but  in  the  opposite  direction 
(Figure  3).  This  is  a  much  more  difficult  situation  since  the  physical  agent  must  move  away 
from  the  goal  to  find  its  best  route. 
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Figure  3  :  Single  Right  Ridge 

4.     Double  Ridge 

This  terrain  has  density  areas  producing  a  right  ridge  on  top  of  a  left  ridge.  Ridge 
densities  are  15  with  a  valley  of  2  in  between  (Figure  4).  An  s-shaped  path  to  get  from  start 
to  goal  is  required  to  transit  this  terrain. 


2   2   2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2   2   2 

2   2   2 

2 

4 

6 

6 

6 

6 

6 

6 

6 

6 

6   6   6 

2   2  4 

6 

8 

a 

a 

a 

a 

a 

a 

a 

a 

a  a  a 

6   6   8 

a 

c 

£ 

f 

f 

£ 

f 

f 

f 

£ 

£   £  £ 

4   6   6 

8 

a 

c 

c 

c 

c 

c 

c 

c 

c 

c   c   c 

4  4   6 

6 

8 

a 

a 

a 

a 

a 

a 

a 

a 

a  a  a 

4  4  4 

6 

6 

8 

8 

8 

8 

8 

8 

8 

8 

8   8   8 

2   2   2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2   2   2 

8   8   8 

8 

8 

8 

8 

8 

8 

8 

8 

6 

6 

4  4  4 

a  a  a 

a 

a 

a 

a 

a 

a 

a 

a 

8 

6 

6  4  4 

c   c   c 

c 

c 

c 

c 

c 

c 

c 

c 

a 

8 

6   6  4 

f   f   f 

f 

f 

f 

f 

f 

£ 

£ 

£ 

c 

a 

8   6   6 

a  a  a 

a 

a 

a 

a 

a 

a 

a 

a 

8 

6 

4   2   2 

6   6   6 

6 

6 

6 

6 

6 

6 

6 

6 

4 

2 

2  2   2 

2   2   2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2   2   2 

2   2  2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2  2   2 

Figure  4  :  Double  Ridge  Terrain 


5.     Single  Left  Plateau 

This  terrain  is  characterized  by  a  large  area  of  high  density  (10)  (dense  but 
passable)  starting  from  left  center  moving  horizontally  out  past  the  grids  midpoint  (Figure 
5).  The  start/goal  density  blocks  are  (4,0)/(l  1,15).  A  successful  transit  can  consist  of  either 
direct  passage  through  the  plateau  or  circumnavigate. 
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Figure  5  :  Single  Left  Plateau  Terrain 


6.     Single  Left  Plateau  with  Ridges 

The  same  as  single  left  plateau,  except  it  has  a  ridge  of  higher  density  (12)  (hard 
to  pass)  along  the  plateaus  perimeter  (Figure  6).  The  start/goal  density  areas  are  (4,0)/ 
(11,15).  Circumnavigation  will  usually  be  the  only  viable  option. 
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Figure  6  :  Single  Left  Plateau  with  Ridges  Terrain 


7.     Random  Terrains 

Four  different  terrains  were  generated  with  random  block  densities  set  anywhere 
from  0  to  15  with  equal  probability.  Shown  in  Appendix  B  (Figures  30,  31,  32  &  33),  these 
grids  were  used  to  test  the  effectiveness  of  the  different  search  strategies  through  unnatural 
terrains. 
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in.  BACKGROUND  WITH  DISCUSSION  PERTINENT  TO  MODEL 

A.  KNOWN  TERRAIN  SEARCH 

A  guarantee  of  optimal  path  from  start  to  goal  is  usually  the  main  concern  in  known 
terrain  search. 

1.     A*  Search 

[Ha68]  -  Widely  accepted  as  the  best  algorithm  for  finding  the  shortest  path  in  a 
known  search  space,  it  uses  the  actual  distance  from  start  and  a  heuristic  estimating  the 
distance  from  goal.  The  object  of  any  frontier  node  is  given  by  the  following  equation: 

/■(«)   =  s(n)+h(n)  (Eq  3.1) 

where  s  (n)  =  the  actual  distance  from  start  to  n,  the  frontier  node,  and  h(n)  =  the 
heuristic  estimated  distance  from  n  to  goal.  Although  guaranteed  to  find  the  optimal  path, 
assuming  that  the  heuristic  estimate  is  always  less  than  the  actual,  it  is  not  required  to 
minimize  the  number  of  nodes  examined.  Using  Euclidean  distance  as  the  distance  to  goal 
heuristic,  A*  search  provided  the  shortest  path  for  each  terrain  used  in  our  experiment  to 
compare  the  effectiveness  of  each  search  scheme. 

B.  UNKNOWN  TERRAIN  SEARCH  BY  A  PHYSICAL  AGENT 

Assuming  limited  sensory  range,  the  physical  agent  cannot  find  the  shortest  path 
without  excessive  moves.  Although  the  shortest  path  would  be  nice,  more  important  to  the 
search  schemes  success  is  the  energy  expended/time  spent  finding  a  satisfactory  path. 
Distance  traveled  is  of  major  concern  as  is  time  to  determine  next  move  (related  to  this  is 
the  computational  limitations  of  the  physical  agent).  The  following  algorithms  have  been 
considered  in  an  attempt  to  find  the  best.  In  all  equations,  n  represents  one  of  the  frontier 
nodes  on  the  frontier  list  unless  otherwise  stated. 


1.  Hill-climb  Search 

[Wi92]  -  Depth  first  search  with  each  move  determined  by  the  best  (least  distance 
from  goal)  adjacent  frontier  node;  n  In  the  following  equation: 

f(n)  =  h(n)  (Eq  3.2) 

When  no  frontier  nodes  are  adjacent  to  the  current,  the  algorithm  backtracks  until  adjacent 
frontier  nodes  are  encountered.  This  search  scheme  minimizes  jumps  in  search  of  the  best 
frontier  node,  but  pays  the  price  in  extra  steps  due  to  unchecked  persistence  on  initially 
good  but  eventually  poor  tracks  and  the  often  resultant  steps  required  to  backtrack. 

2.  Real-time-A*  Search 

[Ko90]  -  Uses  distance  from  current  (actual)  and  distance  from  goal  (heuristic)  to 
determine  best  next  move.  This  search  only  looks  at  adjacent  nodes  (frontier  and  visited). 
In  the  following  equation  n  stands  for  adjacent  non  obstacle  nodes. 

f(n)   =  g(n)+h(n)  (Eq  3.3) 

g  (n)  is  the  actual  distance  from  the  current  to  the  adjacent  node  n.  The  h  (n)  is  a  heuristic 
predicting  the  distance  from  n  to  the  goal.  Initially,  h  (n)  is  calculated  by  using  Euclidean 
distance  in  our  example.  The  algorithm  picks  the  adjacent  node  with  the  best  f(n) .  Before 
moving,  the  value  of  the  f(n)  of  the  second  best  adjacent  node  is  stored  in  the  current  node. 
This  stored  value  will,  in  future  computations,  be  the  node's  h  (n) .  This  value  remains 
constant  until  the  node  is  visited  again.  This  well  conceived  search  scheme  requires 
minimal  computations  and  memory,  but  is  subject  to  wasted  moves  when  drawn  into  local 
traps. 

3.  Best-first  Search 

[Wi92]  (modified  for  physical  agent)  -  Uses  only  the  distance  from  goal  heuristic 
(Euclidean  distance)  to  select  the  next  move. 

f(n)   =  h(n)  (Eq  3.4) 
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Once  selected,  it  uses  the  shortest  path  through  visited  nodes  (known  search  space)  to 
travel  the  distance  from  the  current  node  to  the  selected  frontier  node.  Although,  after  each 
move  it  is  at  the  best  known  location,  the  cost  of  getting  there  can  be  expensive.  In  worst 
case  situations  it  can  end  up  jumping  large  distances  back  and  forth  while  zeroing  in  on  the 
goal. 

4.     Persistence  Search 

[Sh91]  -  Similar  in  concept  to  Real-time-A*,  it  uses  the  distance  from  goal 
heuristic  and  a  weighted  distance  from  current  to  determine  its  next  move.  Unlike  Real- 
time-A*,  it  makes  more  use  of  known  information;  resulting  in  better  moves,  but  decreased 
computational  efficiency.  The  distance  from  goal  is  Euclidean.  The  distance  from  current 
to  frontier  is  the  shortest  path  through  visited  nodes  as  in  Best-first,  but  this  distance  is 
weighted  and  used  in  determining  the  next  move.  The  object  of  a  frontier  node  n  is  given  by 

/(*)  =pfxg(n)+h(n)  (Eq  3.5) 

where  g  (n)  =  shortest  distance  from  current  position  to  n  through  visited  nodes,  h(n)  = 
Euclidean  distance  from  n  to  goal.  A  persistence  factor  (pf  =  0.0  to  1.0)  is  added  to  vary 
the  relative  contribution  of  each  of  the  heuristics  to  the  determination  of  next  move. 
Distance  from  current,  assumed  to  be  always  less  pertinent,  can  be  reduced  in  importance 
in  comparison  to  distance  from  goal. 

C.     GENETIC  ALGORITHMS 

Genetic  algorithms,  developed  by  John  Holland  [Go89]  and  his  associates,  are  based 
on  the  laws  of  natural  selection  and  survival  of  the  fittest.  Subjecting  a  population  (animals, 
search  schemes,  etc.)  to  environments  where  fitness  for  survival  is  required,  individuals 
best  suited  for  survival  will  flourish  and  reproduce  while  individuals  lacking  the  diversity 
required  to  continue  in  all  possible  environments  will  discontinue. 

The  key  to  the  success  of  a  population  is  its  robustness  [Go89].  An  individual,  and 
therefore  a  population,  is  made  up  of  traits  which  are  derived  from  specific  genes  in  the 
individuals  chromosome  [St77].  Applicable  traits  in  the  animal  world  are  weight,  height, 
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leg  length,  neck  length,  etc.  A  combination  of  these  traits  describe  an  individual.  Extremes 
in  any  one  trait  usually  means  more  specification  and  added  survivability  in  a  limited  range 
of  environments,  whereas  moderation  in  traits  means  added  adaptability  for  diverse 
environments.  The  key  is  to  find  the  balance  of  these  two  in  a  population  to  give  it  proper 
robustness.  Example:  the  giraffe  can  afford  to  be  specifically  designed  for  reaching  (long 
neck  and  legs)  because  it  doesn't  face  diversities  in  environment  that  would  require  escape 
through  low  canopy  jungles.  It  is  perfectly  adapted  for  life  on  the  plains  with  occasional 
trees. 

Similarly,  search  strategies  can  be  very  specialized  in  simple  environments.  Search 
through  a  low  density  (of  obstacles)  terrain  can  be  successfully  accomplished  with 
efficiency  and  consistency  by  a  simple  hill  climb  algorithm  (only  one  trait,  distance  to  goal 
of  adjacent  nodes,  is  important).  Search  problems  involving  more  complicated  and 
diversified  solutions  require  the  proper  balance  of  traits  (heuristics)  to  solve.  Simple  direct 
"hill  climbing"  approaches  can  often  result  in  searching  locally  optimal  blind  alleys. 

One  possible  means  of  developing  the  balance  of  traits  required  to  avoid  getting  stuck 
in  the  local  minimum  is  to  enumerate  all  possible  combinations.  This  would  most  assuredly 
find  the  optimum,  but  in  many  problems  the  combinatorial  explosion  of  possibilities  make 
this  method  prohibitive.  Purely  random  combination  of  trials  is  a  possibility  that  seems  to 
avoid  both  the  local  minima  and  the  combinatorial  explosion  problems.  But  on  further 
examination,  it  suffers  the  same  drawbacks  as  enumeration,  in  that  there  are  only  a  limited 
number  of  trials  possible  whether  you  look  at  them  in  order  or  at  random.  Genetic 
algorithms  use  randomness  as  a  tool  in  a  direct  search  for  the  optima.  Promising  potential 
solutions  can  be  searched  in  parallel  while  feedback  information  is  used  to  select  the  next 
partially  random  strategy.  The  results,  as  evaluated  by  DeJong  [De75],  show  the  superiority 
of  genetic  algorithms  over  purely  random. 

The  basic  genetic  algorithm  makes  use  of  a  population  of  individuals  (usually  binary 
strings  of  fixed  length)  that  are  made  up  of  the  traits  pertinent  to  the  problem  (traits  are 
usually  represented  by  a  fixed  number  of  the  bits  in  specific  locations).  Three  genetic 
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operators  are  used  to  transform  the  original  (randomly  generated)  population  into  an 
optimal  one:  selection,  crossover,  and  mutation. 

The  Fitness  of  an  individual  of  the  population  is  established  by  some  form  of 
evaluation  function.  One  scheme  is  to  compare  each  to  a  known  optimum,  assigning  higher 
fitness  to  ones  approaching  the  optima.  This  evaluation  can  also  be  averaged  over  some  set 
number  of  trials  (cycles)  for  each  individual  and  then  assigned  as  the  fitness  before  forming 
the  next  generation. 

Each  new  generation  of  the  population  is  formed  by  stochastically  selecting 
individuals  from  the  prior  population.  Higher  fitness  individuals  have  a  higher  chance  of 
being  selected.  Reproduction  is  performed  by  randomly  pairing  selected  individuals  for 
crossover  and  mutation. 

Crossover  is  performed  at  a  random  point  in  the  binary  string.  The  two  selected  strings 
interchange  their  tail  sections  at  the  crossover  point  to  from  two  new  individuals.  The 
crossover  point  can  be  anywhere  from  0  to  the  last  bit.  For  example,  let  the  two  selected 
strings  be  00000000  and  11111111,  and  let  5  be  the  crossover  point.  Then  the  crossover 
operation  will  produce  the  new  strings  00000111  and  11111000.  In  general,  crossover 
forms  two  new  individuals  with  one  hopefully  having  all  the  best  from  its  two  parents. 

Mutation  is  a  bit  by  bit  operator  that  takes  each  individual  and  randomly  (with  a 
specified  probability)  decides  if  each  bit  will  be  changed?  For  example,  changing  the 
second  bit  of  the  string  000001 1 1  by  mutation  will  result  in  the  new  string  010001 11.  This 
genetic  operator,  as  in  nature,  ensures  that  populations  maintain  adaptability  even  when 
specialization  is  the  rule.  An  extremely  high  mutation  probability  regresses  the  genetic 
algorithms  to  a  uniform  randomly  distributed  population,  a  very  low  one  reduces  the 
populations  adaptability.  A  happy  medium  seems  to  be  in  the  0.01  to  0.001  range  for 
probability  of  individual  bit  mutation. 

A  myriad  of  variations  are  possible  to  improve  the  performance  and  robustness  of  the 
genetic  algorithms.  For  the  purpose  of  our  research,  these  were  found  to  be  unnecessary, 
and  will  not  be  covered  in  this  discussion. 
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IV.  FACTORS  RELEVANT  TO  SEARCH 

A.  DISTANCE  FROM  START 

This  is  usually  the  actual  shortest  path  from  the  start  node  to  the  considered  frontier. 
Currently  believed  to  be  useless  in  a  real-time  environment,  it  should  be  selectively 
eliminated  by  natural  selection  as  the  genetic  algorithm  trains.  For  our  implementation,  it 
is  approximated  by  computing  the  Euclidean  distance  from  start  to  frontier.  It  may  be 
significant  in  some  of  the  more  complex  terrains  that  require  a  switch  back. 

B.  DISTANCE  FROM  CURRENT 

The  distance  from  the  current  node  to  the  frontier  node;  important  in  Real-Time-A* 
and  Persistence  Search  to  determine  if  backtracking  is  worth  the  steps  required.  It  is  the 
actual  distance  computed  as  the  actual  steps  required  to  move  from  the  current  node  to  the 
frontier. 

C.  DISTANCE  FROM  GOAL 

The  Euclidean  distance  from  the  current  node  to  the  goal  node.  This  heuristic  is  usually 
considered  important  in  any  search.  It  is  used  in  combination  with  "distance  from  current" 
for  Persistence  Search,  and  by  itself  for  Best-first  Search. 

D.  CROWDING 

The  crowding  parameters,  crowding  sides  and  crowding  diagonals,  are  an  attempt  to 
assist  the  physical  object  in  avoiding  areas  of  increased  obstacle  density.  This  reduces 
exploration  of  paths  through  high  density  areas,  favoring  the  safer  path  of  increased  options 
available  in  the  open  space.  The  parameters  are  separated  in  case  one  is  more  appropriate 
than  the  other.  Both  would  be  much  more  effective  without  the  self  imposed  constraint  of 
physical  object  perception  only  being  adjacent  nodes.  If  all  nodes  adjacent  to  the  frontier 
node  could  be  seen,  these  factors  importance  would  increase  significantly. 
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1.  Crowding  sides 

This  heuristic  examines  the  frontier  node's  known  horizontal/vertical  neighbors 
to  count  the  number  of  obstacles.  Nodes  with  more  known  obstacle  neighbors  are  less 
desirable.  The  minimum  value  is  0  and  4  is  the  maximum. 

2.  Crowding  diagonals 

This  is  similar  to  the  previous  parameter  with  the  count  being  made  of  the  frontier 
node's  diagonal  neighbors. 

E.  MOVE  AWAY  FACTOR 

It  attempts  to  continually  reduce  the  search  space  by  reducing  desirability  of  nodes 
that  increase  the  x  and/or  y  difference  between  the  current  and  goal  nodes.  Increasing  the  x 
or  y  distance  counts  as  2,  increasing  both  counts  as  4,  and  no  increase  results  in  the  heuristic 
having  a  value  of  0. 

F.  MOMENTUM 

This  heuristic  attempts  to  avoid  zigzag  by  making  forward  (in  relation  to  last  move) 
nodes  the  most  desirable.  It  should  be  useful  in  valley/ridge  terrains  where  the  best  path  is 
straight  through  the  valley.  By  maintaining  momentum,  the  physical  object  avoids  steps 
wasted  in  popping  in  and  out  of  each  crevice  which  has  nodes  closer  to  the  goal.  Straight 
ahead  movement  results  in  a  value  of  0,  a  45°  shift  makes  it  1,  a  90°  shift  is  2,  and  a  135° 
shift  or  non-adjacent  move  results  in  a  value  of  3  (making  only  the  adjacent  nodes  subject 
to  change  after  a  move). 
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V.  PROGRAM  DEVELOPMENT 

A.     DATA  STRUCTURES 

1.  Node  structure 

The  64x64  grid  is  internally  represented  as  a  66x66  two  dimensional  array  (the 
perimeter  nodes  are  all  marked  as  obstacles)  made  up  of  pointers  to  node  records.  The 
records  store  information  pertinent  to  terrain,  search  (heuristics),  graphic  display,  and 
pointers  to  other  node  records  (used  in  the  program  for  various  dynamic  structures).  The 
heuristic  values  stored  include  distance  from  start,  distance  from  goal,  distance  from 
current,  crowding  sides,  crowding  diagonals  and  subtotal.  No  other  node  records  are  used 
in  the  program;  other  structures  requiring  nodes  are  set  up  using  pointers  to  these  records. 

2.  Population  structure 

A  32  member  array  of  individual  records  makes  up  the  population.  Each  stores 
the  individual's  fitness  and  its  chromosome  which  contains  biases  for  each  search 
parameter.  The  chromosome  is  a  32  bit  unsigned  integer;  subdivided  into  eight  four-bit 
unsigned  integers,  it  holds  up  to  eight  heuristic  bias  factors  with  a  range  from  0  to  15. 
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Figure  7  :  Individual  Chromosome  Structure 


16 


3.  Frontier  Heap 

Implemented  as  an  array  of  pointers  to  node  records.  Functions  to  manipulate  this 
min-heap  are  in  theap.c. 

4.  Mate  Heap 

An  array  of  31  integers.  The  first  is  a  count  of  the  heaps  members,  the  other  30 
are  chosen  at  random  in  the  range  of  0  to  the  sum  of  all  individual  fitnesses.  They  are  used 
to  randomly  select  individuals  for  reproduction.  Functions  to  manipulate  this  heap  are  in 
eheap.c. 

B.       SEARCH  ALGORITHMS 

Algorithms  used  for  this  analysis,  with  appropriate  modifications,  are  covered  in  this 
section.  Each  node,  represented  as  a  record,  has  a  number  of  fields  that  are  used  to  store 
needed  information.  When  a  function  operating  on  a  specific  node  v  is  used,  a  read  or  write 
to  the  appropriate  field  takes  place.  For  example:  in  the  A*  algorithm  the  following 
calculation  takes  place;  f(v)  =  s(v)  +  h(v).  The  value  s(v)  calculated  earlier  was  stored  in 
v.s  (the  nodes  dist_from_start  field),  h(v)  is  calculated  and  stored  in  v.h  (the  nodes 
dist_from_goal  field),  and  f(v)  is  stored  in  v.f  (the  nodes  subtotal  field). 

1.  A*  Search 

This  search  (Figure  8)  will  fmd  the  shortest  path  from  start  to  goal  if  a  path  exists. 
The  heuristic  estimating  distance  to  goal  (h(v))  is  a  lower  bound  of  the  actual  cost  of  the 
optimal  path  from  v  to  the  goal. 

2.  Hill-Climb  Search 

This  search  (Figure  9)  always  moves  forward  until  there  is  no  where  to  go.  It  then 
back-tracks  the  way  it  came  until  a  move  is  possible.  It  is  a  depth-first  search  with  a 
heuristic  to  determine  the  best  move  to  advance  it  to  the  goal. 
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a_star_search 

(1) 

current  :=  start 

(2) 

while  current  !=  goal  do 

(3) 

for  all  nodes,  v,  adjacent  to  current  do 

(4) 

if  UNTOUCHED 

(5) 

s(v)  =  currents  +  distance(current,  v)       /*  Euclidean  */ 

(6) 

f(v)  =  s(v)  +  h(v)  /*  h(v)  is  the  Euclidean  dist  to  goal  */ 

(7) 

add  to  frontier  heap 

(8) 

elsif  FRONTIER 

(9) 

if  s(v)  >  currents  +  distance(current,  v) 

(10) 

update  s(v)  and  f(v) 

(11) 

update  position  in  frontier  heap 

(12) 

endif 

(13) 

endif 

(14) 

if  frontier_heap  is  empty 

(15) 

return  BIG_NUMBER  /*  there  is  no  path  from  start  to  goal  */ 

(16) 

endif 

(17) 

end  for  loop 

(18) 

current  :=  top(frontier_heap) 

(19) 

end  while  loop 

(20) 

return  goal.s 

Figure  8  :  A*  Search  Algorithm 
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hillclimbsearch 

(1)  current  :=  start 

(2)  while  current  !=  goal  do 

(3)  best  :=  dummy_node  /*  f(dummy_node)  =  BIG_NUMBER  */ 

(4)  for  all  nodes,  v,  adjacent  to  current  do 

(5)  if  UNTOUCHED 

(6)  f(v)  :=  distance(v,  goal)  /*  the  Euclidean  dist  to  goal  */ 

(7)  mark  v  as  FRONTIER 

(8)  endif 

(9)  if  FRONTIER  and  f(best)  >  f(v) 

(10)  best:=v 

(11)  endif 

(12)  end  for  loop 

(13)  previous_current  :=  current 

(14)  if  best  !=  dummy_node 

(15)  current  :=  best 

(16)  current.predicessor  :=  previous_current 

(17)  elsif  current  !=  start 

(18)  current  :=  current.predicessor 

(19)  else 

(20)  return  BIG_NUMBER 

(21)  endif 

(22)  current.dist_traveled  := 
previous_current.dist_traveled  +  distance(current,  previous_current) 

(23)  end  while  loop 

(24)  return  goal.dist_travekd 

Figure  9  :  Hill  Climb  Search  Algorithm 
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3.  Real-Time-A*  Search 

This  algorithm  (Figure  10)  is  in  accordance  with  Korf' s  description  [Ko90].  For 
our  implementation,  the  node  array  was  used  to  store  the  h  value  since  it  was  already  in 
place,  negating  the  necessity  for  a  hash  table. 

4.  Best-First  Search 

This  search  (Figure  1 1)  always  goes  to  the  best  (minimum  h(v))  node  regardless 
of  its  distance  from  the  current  node.  It  is  possible  to  implement  as  a  specific  case  of  the 
multi-heuristic  search  (Figure  13). 

5.  Persistence  Search 

Shown  in  Figure  12,  gf+hfare  intended  to  effectively  replace/descretize/expand 
the  persistence  factor,  pf,  in  the  original  work  ([Sh91]  equation  3.5).  pf  can  have  any  value 
between  0.0  and  1.0.  We  found  that  an  infinite  range  of  possibilities  for  this  factor  was  not 
required.  A  descrete,  yet  sufficient,  span  can  be  obtained  by  setting  gf  and  hf  to  any 
number  of  possibilities  where  gf<  hf.  Setting  hf  to  15  and  incrementing  gf  from  0  to  15 
gives  us  the  equivalent  of  a  0.0  to  1.0  range  incriminated  by  0.067. 

/(v)   =  gfxg(v)+hfxh(v)  (Eq  5.1) 

There  is  also  now  the  expanded  capability  of  having  the  g(v)  be  the  more 
important  factor  in  the  search  (gf>  hf).  This  search  can  also  be  implemented  as  a  specific 
case  of  the  multi-heuristic  search. 
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real_time_astar_search 

(1)  current  :=  start 

(2)  best  :=  dummy_node 

(3)  second_best  :=  dummy_node  /*  f(dummy_node)  =  BIG_NUMBER  */ 

(4)  while  current  !=  goal  do 

(5)  for  all  nodes,  v,  adjacent  to  current  do 

(6)  if  UNTOUCHED 

(7)  h(v)  :=  distance(v,  goal)  /*  Euclidean  */ 
/*  else  h(v)  is  already  set  */ 

(8)  endif 

(9)  g(v)  :=  distance(current,  v) 

(10)  f(v):=g(v)  +  h(v) 

(11)  ifbest.f>f(v) 

(12)  second_best  :=  best 

(13)  best:=v 

(14)  elsif  second_best.f  >  f(v) 

(15)  second_best  :=  v 

(16)  endif 

(17)  end  for  loop 

(18)  previous_current  :=  current 

(19)  current  :=  best 

(20)  previous_current.h  :=  second_best.f 

(21)  current.dist_traveled  := 
previous_current.dist_traveled  +  distance(current,  previous_current) 

(22)  end  while  loop 

(23)  return  goal.dist_traveled 

Figure  10  :  Real  Time  A*  Search  Algorithm 
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best_first_search 

(1)  current  :=  start 

(2)  while  current  !=  goal  do 

(3)  for  all  nodes  adjacent  to  current  do 

(4)  if  UNTOUCHED 

(5)  h(v)  :=  distance(node,  goal)  /*  Euclidean  */ 

(6)  add  to  frontier_heap 

(7)  endif 

(8)  end  for  loop 

(9)  if  empty(frontier_heap) 

(10)  return  BIG_NUMBER  /*  no  solution  */ 

(11)  endif 

(12)  v.dist_traveled  :=  current.dist_traveled  +  g(v) 

where  g(v)  is  the  shortest  distance  through  known  paths  from 
current  to  frontier  node. 

(13)  previous_current  :=  current 

(14)  current  :=  top(frontier_heap)  /*  minimum  h(v)  */ 

(15)  end  while  loop 

(16)  return  goal.dist_traveled 

Figure  1 1  :  Best  First  Search 
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persistencesearch 

(1)  current  :=  start 

(2)  while  current  !=  goal  do 

(3)  for  all  nodes  adjacent  to  current  do 

(4)  if  UNTOUCHED 

(5)  h(v)  :=  distance(node,  goal)  /*  Euclidean  */ 

(6)  add  to  frontier_heap 

(7)  endif 

(8)  end  for  loop 

(9)  if  empty(frontier_heap) 

(10)  return  BIG_NUMBER     /*  no  solution  */ 

(11)  endif 

(12)  find  the  frontier  node,  v,  that  minimizes  the  equation: 

f(v)  :=  gf  *  g(v)  +  hf  *  h(v)      where  g(v)  is  the  shortest  distance  through 
known  paths  from  current  to  frontier  node,  gf  and  hf,  set  before  search, 
are  bias  factors  used  to  vary  the  relative  importance  of  g(v)  and  h(v). 
They  can  have  a  value  from  0  to  15. 

( 1 3)  v.dist_traveled  :=  current.dist_traveled  +  g(v) 

(14)  current:=v 

(15)  remove  current  from  frontier_heap  and  update 

(16)  end  while  loop 

(17)  return  goal.dist_traveled 

Figure  12  :  Persistence  Search  Algorithm 
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6.     Multiheuristic  Search 

This  is  the  general  algorithm  (Figure  13)  enstantiated  in  our  case  to  handle  five 
stablejieuristics  and  two  unstable_heuristics.  Stablejieuristics  being  ones  that  have 
values  that  will  not  change  if  more  than  two  steps  away  from  the  current  node.  They  include 
Euclidean  distance  from  goal  (hg),  Euclidean  distance  from  start  (hs)>  crowd  sides  (hcs), 
crowd  diagonals  (hcd),  and  momentum  (hm).  The  subtotal  fs  (v)  is  calculated  using  these 
functions  multiplied  by  their  respective  bias  factor  and  stored  in  v. subtotal. 
fs(v)  =  hgfxhg(v)  +  hsfxhs(v)  +  hcsfxhcs(v)  +  hcdfxhcd(v)  +  hmfxhm(v)    (Eq  5.2) 

Unstable_heuristics  have  values  that  are  liable  to  change  as  the  current  node 
changes.  Examples  in  our  case:  distance  from  current  (hdc)  and  move  away  (hma).  The 
algorithm  minimizes  equation  5.3  using  the  efficient  "branch-and- bound"  search  through 
known  (visited)  nodes  described  in  section  4.3  of  [Sh91]. 

/(v)  =  fs  (v)  +  hdcfx  hdc  (v)  +  hmafx  hma  (v)  (Eq  5.3) 

The  hsf,  hgf,  hdcf,  hcsf,  hcdf,  hmaf,  and  hmf  are  bias  factors  that  correspond  with  the 
individual  chromosome's  lower  28  bits  which  are  set  during  training.  The  highest  four  bits 
are,  in  our  implementation,  a  place  holder  for  future  additional  heuristics  since  only  seven 
applicable  heuristics  were  identified.  Note  that  the  Best-first  and  Persistence  Search  could 
be  implemented  as  special  cases  of  the  multi-heuristic  search  algorithm.  Best-first  uses  an 
individual  chromosome  input  of  00100000  (the  third  factor  being  hgf).  Persistence  Search 
uses  an  individual  chromosome  input  of  OOxyOOOO  with  x  and  y  varying  from  0  to  15 
(fourth  factor  being  hdcf). 
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multi_heuristic_search 

(1)  current  :=  start 

(2)  while  current  !=  goal  do 

(3)  for  all  nodes  v  within  2  moves  of  current  do 

(4)  if  adjacent  and  UNTOUCHED 

(5)  v.subtotal  :=  inner_product(stable_heuristics  *  respective_biases) 

(6)  add  v  to  frontier_heap         /*  min  subtotal  node  on  top  */ 

(7)  elsif  FRONTIER 

(8)  if  any  stable_heuristics  of  v  have  changed 

(9)  v.subtotal  :=  v.subtotal  +  adjustment 

(10)  update  position  in  frontier_heap 

(11)  endif 

(12)  endif 

(13)  end  for  loop 

(14)  if  empty  (frontier_heap) 

(15)  return  BIG_NUMBER     /*  no  solution  */ 

(16)  endif 

(17)  find  frontier  node,  v,  that  minimizes 

f(v)  =  v.subtotal  +  inner_product(unstable_heuristics  *  respective_biases) 

( 1 8)  v.dist_traveled  :=  current.dist_traveled  +  g(v) 

where  g(v)  is  the  shortest  distance  through  known  paths  from  current  to 
frontier  node. 

(19)  current  :=  v         /*  and  remove  v  from  heap  */ 

(20)  end  while  loop 

(21)  return  goal.dist_traveled 

Figure  13  :  Multi  Heuristic  Search  Algorithm 
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C.     GENETIC  ALGORITHM 

The  task  of  the  genetic  algorithm  is  to  find  the  combination  of  the  seven  bias  factors 
that  will  result  in  the  optimum  search  scheme.  The  values  of  these  seven  bias  factors  are 
stored  in  a  single  individuals  chromosome.  Application  of  genetic  operators  to  a  population 
(32  in  our  case)  of  these  individuals  will,  after  numerous  iterations,  produce  our  desired 
optimal  individual. 

The  genetic  algorithm,  described  in  this  section,  is  invoked  during  training  after  some 
predetermined  number  of  cycles  (making  up  one  generation).  The  input  population  will 
have  a  fitness  value  (ability  to  get  through  the  terrain)  assigned  to  each  of  it's  32  individuals 
(details  of  this  process  are  described  in  the  next  chapter).  This  fitness  value  and  the 
individual's  chromosomal  make-up  are  required  by  the  genetic  algorithm. 

Our  algorithm  (Figure  16)  makes  use  of  the  three  genetic  operators:  selection, 
crossover,  and  mutation.  The  implementation  is  similar  to  the  algorithm  presented  in 
chapter  one  of  the  text  by  Goldberg,  [Go89],  with  the  additions  of  allowing  the  best  two 
individuals  to  go  unchanged  and  an  average  of  one  out  of  seven  of  the  remaining  not  going 
through  crossover.  The  result  is  similar  to  De  Jong's  R3  elitist  model  [De75].  Examples  of 
our  crossover  implementation  are  detailed  in  Figures  14  &  15. 


Alleles  are  represented  in  hexidecimal 

Before  allele  crossover:  55555555  /  88888888 

Randomly  picked  crossover  allele  position  is  3  (4th  allele) 
After  allele  crossover:  55558888  /  88885555 


Figure  14  :  Allele  Crossover  Example 


The  4th  allele  is  expanded  out  into  binary  representation 

Before  bit  crossover:  555  0101  8888  /  888  1000  5555 

Randomly  picked  crossover  position  between  bits  is  2 

After  bit  crossover:  555  1001  8888  /  888  0100  5555 


Figure  15  :  Bit  Crossover  Example 
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geneticalgorithm 

(input  is  a  population  of  individuals) 

(1)  total_fitness  :=  all  individual  fitnesses  added  together 

(2)  select  32  individuals  as  follows  /*  selection  */ 

(3)  best  :=  individual  with  the  highest  fitness 

(4)  second_best  :=  individual  with  the  second  highest  fitness 
/*second_best  must  be  distinct  from  best  */ 

(5)  stochastically  select  30  individuals  with  higher  fitness  individuals  having 
the  greatest  chance  of  selection 

(6)  end  selection 

(7)  create  new_population  with  these  32  individuals 

pair  individuals  in  such  a  way  that  it  is  unlikely  that  an  individual  is  paired 
with  itself;  pair  best  with  second_best 

(8)  for  each  individual  pair,  except  best  and  second  best,  do 

(9)  randomly  pick  crossover  allele  position  /*  crossover  */ 

(10)  if  not  0     /*  0  means  no  crossover  */ 

(11)  exchange  all  alleles  after  the  crossover  allele 

(12)  randomly  pick  crossover  position  between  bits  of  selected  allele 

(13)  if  not  0  or  4     /*  0  or  4  means  crossover  does  not  breakup  the  allele  */ 

(14)  exchange  bits  after  the  crossover  position  between  bits 

(15)  endif 

(16)  endif 

(17)  for  each  gene  of  the  individual  /*  mutation  */ 

(18)  invert  bit  if  random  <  prob  of  mutate 

(19)  end  for  loop 

(20)  end  for  loop 

(21)  add  best  +  second_best  to  new_population  as  individuals  0  &  1  respectively. 

(22)  return  new_population 

Figure  16  :  Genetic  Algorithm 
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VI.  TRAINING 

Training  of  the  population  is  analogous  to  selectively  breeding  a  random  group  of 
asexual  organisms  to  obtain  superior  capability  in  a  specific  area.  The  capability  you  wish 
to  optimize  is  transit  from  start  to  goal  in  the  least  number  of  steps.  The  specific  area  is  a 
specific  terrain  layout  where  you  have  an  idea  about  general  areas  of  density,  but  have  no 
information  about  the  location  of  specific  obstacles. 

The  first  step  is  to  generate  a  series  of  specific  terrains  from  your  general  idea  of  the 
densities.  This  can  be  done  by  placing  obstacles  in  each  area  if  a  randomly  generated 
number  is  less  than  the  specified  density.  In  our  implementation,  we  simply  loop  through 
the  64x64  node  array  assigning  each  nodes  state  to  OBSTACLE  if  the  random  number  is 
less  than  the  density  value  of  the  corresponding  density  block.  The  second  step  is  to 
generate  a  population  of  32  individuals  giving  them  randomly  generated  chromosomes. 
Now  the  training  begins  (Figure  17).  In  all  our  work,  we  used  1000  generations  with  five 
cycles  (trials)  per  generation. 

The  returned  population  evolves  through  the  trials  of  5000  terrains.  One  of  the 
individuals  of  this  population  is  likely  to  have  a  chromosome  that  approximates  the 
optimum  combination  of  bias  factors.  Identification  of  this  individual  is  accomplished 
during  testing. 
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training 

( 1 )  for  the  number  of  generations  do 

(2)  for  the  number  of  cycles  do 

(3)  loop  until  a  successful  A*  search 

(4)  create  a  terrain  from  the  density _array 

(5)  shortest_path  :=  A*  search 

(6)  end  until  loop 

(7)  run  each  individual  through  the  terrain  accumulating  its  fitness_sum  by 
comparing  its  actual  path  to  the  shortest  path 

(8)  end  for  loop 

(9)  compute  each  individual's  average  fitness  from  fitness_sum  and 
number  of  cycles 

(10)  apply  the  genetic  algorithm  to  the  population 

(11)  end  for  loop 

(12)  return  a  trained  population. 

Figure  17  :  Training  Algorithm 
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VII.  TESTING 

Testing  of  the  trained  populations  was  performed  by  comparing  the  search  conducted 
by  the  best  individual  in  each  population  to  searches  accomplished  using  Hill-climbing, 
Best-first,  Real-time- A*,  and  Persistence  search.  The  following  equation  was  used  to 
compute  fitness  for  all  search  schemes: 

fitness  =  integer (( (shortestpath)  +  (actualpath))  x  100)  (Eq  7.1) 

Each  search  scheme  was  tested  on  500  distinct  terrains  produced  using  the  corresponding 
density  matrix. 

Before  testing,  the  best  of  each  population  was  chosen  by  running  the  population 
through  50  distinct  terrains.  The  individual  with  the  highest  fitness  was  chosen  to  represent 
the  GA-trained  population.  The  best  values  for  distance  from  goal  and  distance  from 
current  bias  factors  for  the  Persistence  search  were  determined  by  running  32  combinations 
(chromosomes  of  OOfOOOOO  to  OOffOOOO  and  OOOfOOOO  to  OOffOOOO)  through  50  distinct 
terrains.  From  this,  the  best  combinations  of  the  two  factors  was  used  to  represent 
Persistence  search. 

The  GA-produced  best  individual,  Persistence  best,  Hill-climb,  Real-time  A*,  and 
Best-first  schemes  were  then  all  used  to  find  a  path  in  the  500  separate  terrains.  Average 
fitnesses  over  the  500  were  assigned  and  a  comparison  of  these  fitnesses  is  presented  in  the 
results. 
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VIH.  EXPERIMENTAL  RESULTS 

The  fitness  of  each  search  scheme  in  these  results  is  the  number  of  its  required  steps 
divided  by  the  minimum  steps  possible,  averaged  over  the  500  terrains  used  for  testing. 
Fitness  is  presented  as  a  percentage,  with  a  100%  search  scheme  being  one  that  can,  on  the 
average,  search  a  terrain  type  in  the  minimum  steps  possible.  In  general,  the  easier  the 
density  layout  of  the  terrain,  the  higher  the  fitness  will  be. 

A.      NATURAL  TERRAINS 

A  graph  comparing  the  fitness  of  applicable  search  schemes  is  presented  for  each 
natural  terrain  density  layout  (Figures  18  -  23).  The  following  discussion  is  pertinent  to 
each  of  these  comparisons. 

1.  Central  Mountain 

This  graph  (Figure  18)  shows  that  this  terrain  is  only  moderately  hard  for  all  the 
search  schemes.  Persistence  search  with  a  distance  to  current  factor  (gf)  of  15  and  a  distance 
to  goal  factor  (hf)  of  1 1  (gf/hf  =  15/11)  was  the  best  of  the  conventional  search  methods. 
The  genetic  algorithm  produced  an  individual  with  chromosomal  make-up  of  f00732b9 
(see  figure  7,  page  16  for  breakdown)  which  performed  1.20  times  better  than  the  best 
conventional.  Driven  more  to  the  goal  by  the  move-away  heuristic  than  distance  to  goal, 
this  scheme  was  better  equipped  to  avoid  the  congestion  of  the  central  mass. 

2.  Single  Left  Ridge 

Overall  this  terrain  was  a  little  harder  than  the  Central  Mountain  but  was  still 
handled  moderately  well  by  all  search  strategies  (Figure  19).  The  best  conventional  was 
again  persistence  search  using  a  gf/hf  ratio  of  15/6.  The  genetic  algorithm  scheme 
(f00828ff)  had  a  fitness  1.16  times  as  good  as  the  best  Persistence  and  1.28  times  better  than 
the  next  competitor  (Hill-climbing). 
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Figure  18  :  Central  Mountain  Results 


Best  First 


Persistence 


Hill  Climb 


Real-time-A* 


GA-produced 


FITNESS  (percent) 


Figure  19  :  Single  Left  Ridge  Results 
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3.     Single  Right  Ridge 

As  could  be  predicted,  this  was  a  hard  problem  (Figure  20).  The  normal  search 
schemes  tend  to  spend  a  lot  of  time  searching  every  possible  route  that  was  most  direct  to 
the  goal.  They  would  get  stuck  under  the  ridge  with  no  way  around,  except  back  the  same 
way  they  came.  The  best  of  these  was  Hill-climbing  since  it  probably  doesn't  waste  a  lot  of 
steps  backtracking.  The  best  Persistence  search  had  a  gf/hf  ratio  of  2/15  showing  it's 
favoritism  for  a  no  backtrack  approach.  The  search  scheme  produced  by  the  genetic 
algorithm  was  superior  to  all  by  a  multiplication  factor  of  1.26.  Its  chromosomal  make-up 
was  f00c2ca8.  This  scheme  considers  distance  to  goal  to  be  not  significant.  It  instead  uses 
move  away  factor  as  the  drive  toward  the  goal.  As  the  Persistence  search  with  its  gf/hf  ratio 
of  2/15  the  genetic  algorithm  produced  scheme  considered  the  amount  of  backtracking  a 
major  factor. 
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Figure  20  :  Single  Right  Ridge  Results 
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4.     Double  Ridge 

This  terrain  was  a  very  difficult  problem  (Figure  21).  Requiring  negotiation 
around  two  ridges  which  involved  a  switch-back  away  from  the  goal,  none  of  the  search 
schemes  were  over  20%  fitness.  The  average  of  the  five  schemes  was  10%.  The  best 
Persistence,  with  a  gh/hf  ratio  of  6/15,  was  roughly  equivalent  to  Hill-climbing.  The  genetic 
algorithm  generated  scheme  with  a  chromosomal  make-up  of  f83bl9bc  was  the  best 
strategy  with  a  fitness  1.21  times  better  than  Persistence.  Here  is  an  example  where  distance 
to  start  was  of  significance;  probably  helpful  in  influencing  the  search  to  make  the  switch- 
back away  from  the  goal.  Move  away  factor  was  a  major  influence  in  striving  toward  the 
goal,  backtracking  was  determined  to  be  non-productive,  but  maintaining  momentum  was 
found  to  be  important.  It's  interesting  to  note  that  diagonal  crowding  was  considered  more 
important  than  side  crowding  (no  explanation).  The  complexity  of  this  scheme  with  the 
subtle  interaction  between  these  differing  bias  factors  helps  to  confirm  the  necessity  of  a 
genetic  algorithm  to  sort  them  out. 
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Figure  21  :  Double  Right  Left  Ridge  Results 
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5.     Single  Left  Plateau 

This  terrain  was  slightly  easier  for  all  the  search  strategies  (Figure  22)  although  it 
presented  a  unique  problem.  The  through  the  plateau  route  is  possible  but  requires 
numerous  explorations.  Circumnavigating  the  plateau  saves  exploration  steps  but  costs  in 
the  distance  required.  Since  each  of  the  500  terrains  had  varying  obstacle  placement,  we 
suspect  sometimes  it  was  best  to  transit  through  and  other  times  better  to  go  around.  Since 
no  general  path  was  consistently  optimal,  the  genetic  algorithm  had  to  develope  a  scheme 
that  was  equally  effective  for  both  routes  or  concentrate  on  perfecting  one.  In  either  case, 
its  performance  was  again  superior  by  a  significant  margin  (multiplicative  factor  of  1.19). 
The  resultant  chromosomal  make-up  was  f05e884f.  The  next  best  was  Persistence  with  a 
gf/hf  ratio  of  11/15. 
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Figure  22  :  Single  Left  Plateau  Results 
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6.     Single  Left  Plateau  With  Ridges 

This  terrain  adds  topological  characteristics  that  favor  circumnavigation  as  a 
search  strategy.  The  genetic  algorithm  produced  scheme,  with  a  chromosomal  make-up  of 
f07c033d,  was  the  best  by  a  multiplication  factor  of  1.17  over  the  next  best  competitor 
(Figure  23).  Momentum  being  the  most  important  factor,  it  probably  helped  keep  the  search 
moving  horizontally  until  clear  of  the  plateau.  Persistence  was  again  the  second  best  with 
a  gh/hf  ratio  of  11/15. 
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Figure  23  :  Single  Left  Plateau  With  Ridges  Results 

7.     General  Comment 

The  genetic  algorithm  was  extremely  successful  in  producing  the  best  search 
strategies  for  all  natural  terrains. 
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B.      RANDOM  TERRAINS 

Although  the  genetic  algorithm  produced  search  schemes  where  clearly  superior  for 
the  natural  terrains,  we  wanted  to  test  their  viability  on  randomly  generated  terrains. 

1.     Random  One  /  Random  Two  /  Random  Three 

The  results  from  these  three  terrains  showed  that  the  search  heuristics  produced 
by  genetic  algorithms  was  of  minimal  value  (Figures  24  to  24).  These  were  all  simple 
problems  with  the  average  fitness  for  all  the  search  schemes  being  64%.  Fitness  varied  little 
between  search  strategies  with  a  maximum  of  a  8%  difference  between  the  best  and  the 
worst.  The  genetic  algorithm  produced  scheme  was  1.02  (Random  One),  1.01  (Random 
Two),  and  1 .03  (Random  Three)  times  as  good  as  the  best  conventional  search  strategy.  The 
Random  One  persistence  gf/hf  ratio  was  15/11;  the  genetic  algorithm  produced 
chromosome  was  fle90234.  The  Random  Two  gf/hf  =  15/4;  GA-produced  chromosome  = 
f0b947cl.  The  Random  Three  gh/hf  =  15/5;  GA-produced  chromosome  =  flf73351. 
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Figure  24  :  Random  One  Results 
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Figure  25  :  Random  Two  Results 
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Figure  26  :  Random  Three  Results 
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2.     Random  Four 

This  terrain  was  significantly  harder  than  the  other  three  randomly  generated 
terrains  which  can  be  observed  by  the  low  performance  of  the  search  strategies.  The 
average  fitness  of  all  strategies  was  36%  (Figure  24).  The  difficulty  probably  comes  from 
the  encapsulation  of  the  goal.  Examining  figure  33,  page  46,  we  can  see  that  the  goal  is 
blocked  by  mostly  high  density  blocks  from  (11,15)  down  to  (1 1,10)  across  to  (15,10).  The 
only  passible  blocks  are  (15,10)  and  (11,12).  Neither  of  which  are  a  direct  route, 
necessitating  significant  exploration.  The  genetic  algorithm  produced  (f0b51535)  scheme 
was  1.10  times  better  than  the  next  best  which  was  a  persistence  strategy  with  a  gf/hf  ratio 
of  14/15.  This  again  seems  to  suggest  that  the  genetic  algorithm  is  only  required  when  the 
problem  is  hard. 
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Figure  27  :  Random  Four  Results 
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3.     General  Comments 

It  is  difficult  to  improve  on  the  simple  search  strategies  when  the  terrain  is  of  low 
complexity.  As  suggested  by  DeJong  [De92]  the  genetic  algorithm  can  only  optimize  to  a 
certain  point  (dependent  on  implementation)  before  reaching  a  state  of  dynamic 
equilibrium.  The  first  three  random  terrains  were  of  insufficient  complexity  to  allow  the 
genetic  algorithm  to  convincingly  surpass  all  conventional  search  schemes.  It  was  however, 
in  all  cases,  better  than  the  best  conventional  ones. 

C.     GENERAL  COMMENTS 

In  all  cases,  although  the  genetic  algorithm  produced  strategy  was  always  as  least  as 
good  as  the  next  best,  it  was  not  a  substantial  improvement  over  Persistence  search  unless 
the  terrain  was  natural.  Only  in  the  most  complex  of  the  four  random  terrains  did  the  genetic 
algorithm  produced  scheme  really  excel.  This  seems  to  suggest  that  the  additional 
heuristics  are  only  essential  in  natural  terrains  where  some  pattern  in  obstacle  density  exists 
or  in  random  terrains  of  high  complexity. 

Actual  natural  terrains,  although  usually  best  modeled  by  our  natural  terrains,  could 
possibly  be  more  similar  to  the  random.  Since  the  genetic  algorithm  produced  search 
strategies  are  substantially  better  for  our  natural  terrains  and  as  least  as  good  as  standard 
search  schemes  for  random  terrains,  they  should  be  advantageous  to  use  on  any  actual 
natural  terrain.  This  is  of  course  contingent  on  the  physical  agent's  dependence  on  minimal 
steps  and  its  computational  speed.  If  it's  computational  speed  is  sufficient  to  avoid  delays 
before  each  step  and/or  minimal  steps  are  essential,  the  genetic  algorithm  produced  scheme 
should  always  be  used. 

Appendix  C  shows  a  comparison  of  the  average  time  required  for  each  strategy  to 
search  from  start  to  goal  for  each  of  the  terrains.  As  expected,  the  more  complicated 
strategies  require  additional  computation  time,  but  are  not  considered  slow  enough  to 
prohibit  their  use  except  in  cases  of  high  speed  agents  with  slow  computational  speed. 


40 


IX.  CONCLUSIONS 

Heuristics  previously  used  for  search  of  an  unknown  space  by  a  physical  agent  are 
distance  from  goal  and  distance  from  current.  These  are  insufficient  to  minimize  energy 
expenditure  (steps  taken)  when  some  general  knowledge  of  the  area  is  known.  The 
additional  heuristics  found  to  be  pertinent  are  distance  from  start,  crowding  factors  which 
account  for  obstacle  node  density  around  the  considered  frontier  node,  move-away  factor 
which  encourages  reduction  of  the  search  space,  and  momentum  which  avoids  wasted  steps 
in  course  variations.  These  seven  heuristics  with  their  proper  individual  biases  were  found 
to  be  superior  to  standard  search  schemes.  In  this  thesis  we  showed  that  genetic  algorithms 
can  be  effectively  used  to  develop  optimal  heuristic  biases  that  are  adaptable  to  unknown 
search  spaces  if  some  general  knowledge  of  the  search  space  is  available.  Training  done 
with  randomly  generated  search  spaces  having  common  characteristics  lead  to  robust 
search  schemes  which  are,  on  the  average,  more  fit  than  previously  used  strategies. 

We  believe  that  this  methodology  of  identifying  all  possible  heuristics,  fitting  them 
into  a  binary  representation,  and  applying  genetics-based  training  is  also  applicable  to  a 
multitude  of  real-time  search/optimization  problems.  Tests  in  other  specific  areas  are 
needed  to  prove  our  conjecture.  In  addition,  further  research  could  be  done  in  the 
application  of  more  advanced  genetic  algorithms.  Our  results  showed  significant 
improvement  using  only  basic  genetics-based  concepts,  advanced  techniques  should 
continue  to  improve  the  effectiveness  of  resultant  strategies. 
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APPENDIX  A 
TERRAIN  DEVELOPMENT  FROM  A  SAMPLE  DENSITY  MATRIX 

Terrains  are  randomly  produced  using  a  density  matrix  as  a  guide.  Figure  28  shows  a  density 
matrix  that  was  used  to  develope  the  terrain  is  shown  in  figure  29.  This  density  matrix  was  not  used 
for  our  analysis,  but  helps  to  make  clear  the  relationship  between  the  density  matrix  and  the  actual 
terrain. 

The  density  matrix  is  stored  as  a  text  file  as  shown  in  the  figure.  At  each  cycle  for  training  or 
iteration  for  testing  the  density  information  is  used  to  form  a  new  terrain.  Each  hexadecimal  number 
represents  the  desired  density  for  a  4x4  area.  The  actual  obstacle  placement  is  random.  Compare 
figures  28  and  29.  The  top  left  4x4  area  was  filled  in  by  checking  if  a  random  number  (between  0 
and  15)  is  less  than  4  at  each  node.  This  should  on  the  average  happen  4  out  of  16  times  making  the 
obstacle  count  of  each  4x4  area  equal  4.  The  top  left  4x4  is  the  average  case  with  4  out  of  the  16 
nodes  being  obstacles. 

The  remaining  4x4  areas  are  filled  out  in  a  similar  fashion.  The  start  and  goal  nodes  are  chosen 
at  random  in  the  (2,2)  and  (13,13)  areas  as  also  demonstrated  in  figure  29. 
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Figure  28  :  Sample  Density  Matrix 
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Terrain  has  a  density  of  4  in  the  top  1 1/16  Goal 


Start 


The  lower  5/16  shows  increasing  density  areas  0  to  15 
surrounded  by  areas  of  0  density. 


Figure  29  :  One  of  Many  Possible  Resultant  Terrains 
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APPENDIX  B 
RANDOM  TERRAIN  DENSITY  MATICES 
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Figure  30  :  Random  Terrain  One 
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Figure  31  :  Random  Terrain  Two 
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Figure  32  :  Random  Terrain  Three 
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Figure  33  :  Random  Terrain  Four 
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APPENDIX  C 


SEARCH  STRATEGY  COMPUTATIONAL  TIME 


TABLE  1:    AVERAGE  SECONDS  REQUIRED  TO  SEARCH  EACH  TERRAIN 


Best  First 

Persistence 

Hill  Climb 

Real-Time- 
A* 

GA- 

Produced 

Central 
Mountain 

0.0239 

0.0274 

0.0065 

0.0163 

0.0615 

Single  Left 
Ridge 

0.0232 

0.0264 

0.0075 

0.0176 

0.0694 

Singe  Right 
Ridge 

1.3051 

0.1496 

0.0199 

0.1405 

0.2718 

Double 
Ridge 

3.2553 

1.3562 

0.0420 

0.4131 

2.3420 

Single  Left 
Plateau 

0.1071 

0.0982 

0.0167 

0.0626 

0.1481 

Single  Left 

Plateau 
With  Ridges 

0.0769 

0.0856 

0.0146 

0.0487 

0.1310 

Random 
One 

0.0177 

0.0265 

0.0064 

0.0114 

0.0478 

Random 
Two 

0.0313 

0.0267 

0.0056 

0.0150 

0.0644 

Random 
Three 

0.0294 

0.0258 

0.0065 

0.0151 

0.0521 

Random 
Four 

0.1629 

0.2255 

0.0165 

0.0481 

0.1563 
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APPENDIX  D 
PROGRAM  C  CODE 

I .  ga_search .  h 5C 

2  .  main .  c 5f 

main 

3  .  train  .  c 56 

train 

put_gen 

put_rs 

4.  test.c 61 

test 

5  .  tsetup .  c 65 

get_seed 

get_pers_chrom 

read_density_f ile 

make_array 

make_node 

f ind_node 

6 .  tpopulation .  c 70 

create_population 
ne w_indi vidua 1 
get^population 
put_population 

7.  astar.c 73 

a_star 
update_astar_f rontier 

8.  hill_climb.c 7  6 

hill_climb 
move_ad j  acent 
f ind_best 

9.  rt_astar.c 79 

rta_star 

updat e_rt as tar_adj acent 

insert 

10  .  bfsearch.  c 83: 

bf search 

bf_update_f rontier_list 

bf_pick_best_f rontier 

II .  psearch.  c 87 

psearch 

p_update_f rontier_list 

p_pick_best_f rontier 

12.  tsearch.c 91 

search 
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13.  tfrontier.c 94 

update_f rontier_list 

update_adjacent_obstacles 

pick_best_f rontier 

reset_back_track_state 

update_list 

diff_int 

adjacent 

upda  t  e_c  r  o  wd_s  i  de  s 

update_crowd_diag 

c  a 1 c_mo ve_a wa y 

c a 1 c_mome n t um 

compute_subtotal 

14  .  theap.c 102 

insert_heap 
delete_heap 
move_heap 
swap 

15.  evolve. c 104 

evolve 

create_mate_heap 

allele_crossover 

bit_crossover 

crossover 

get_mask 

mutate 

one_if_mutate 

get_odd 

set-equal 

16.  eheap.c 112 

insert_mate_heap 
pop_mate_heap 
mo  ve_ma  t  e_h  e  a  p 
swap_num 

17  .  tmisc.  c 114 

gen_xi 
gen^yi 
equalf 
show_least_nodes 

18.  tdisplay.c 117 

initialize 

draw_terrain 

show_mouse 

draw_nodes 

draw_grid 

squaref 

square 

19.  tprint.c 122 

print_density 
print_population 
print  node 
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ga_search.h 

/* 

File:  ga_search.h 

Programmer:  g.b.  parker 

Environment:  any 

Language :  C 

Date:  9  July  92 

Revised: 

Comments:  This  file  holds  all  the  header  information  needed  for  all 
ga_search  files. 
*/ 

♦include  <stdio.h> 


#def ine 

F  0 

#def ine 

T  1 

♦define 

UNTOUCHED 

0 

♦define 

OBSTACLE 

1 

♦define 

VISITED 

2 

♦define 

FRONTIER 

3 

♦define 

START 

4 

♦define 

GOAL 

5 

♦define 

CURRENT 

6 

♦define 

SHORTEST 

7 

♦define 

X 

8 

♦define  NUM  0 


♦define 

N 

1 

♦define 

E 

2 

♦define 

S 

3 

♦define 

W 

4 

♦define 

NE 

1 

♦define 

SE 

2 

♦define 

SW 

3 

♦define 

NW 

4 

♦define  MASKO  0  /*  0000  */ 

♦define  MASK1  8  /*  1000  */ 

♦define  MASK2  12  /*  1100  */ 

♦define  MASK3  14  /*  1110  */ 

♦define  MASK4  15  /*  1111  */ 
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♦define  NA  -1 

dex  */ 

#define  BIG_NUMBER        10000 
♦define  STANDARD_DENSITY   4 
♦define  PROB  BIT  MUTATE   50 


/*  not  applicable;  for  frontier  in- 


/*   x/10000  prob  of  mutate  */ 


♦define  SQRT2 


1.414213562 


♦define  randl6()  ( (random( ) /13)  %  16)  /*  return  rand  int  from  0  to  15  */ 

♦define  rand8()   ( (random( ) /13)  %  8)  /*  return  rand  int  from  0  to  7  */ 

♦define  rand5()   ( (random( ) /13)  %  5)  /*  return  rand  int  from  0  to  4  */ 

♦define  randlOOOOO  (random()  %  10000)  /*  return  rand  int  from  0  to  9999  */ 

/*  Graphics  definitions  */ 
♦define  SLEEPTIME    2 


♦define  ASTITLE 
♦define  GATITLE 
♦define  PERSTITLE 
♦define  BFTITLE 
♦define  HCTITLE 
♦define  RTASTITLE 


"A*  Search 

"Genetic  Algorithm  Produced  Search 

"Persistence  Search 

"Best  First  Search 

"Hill  Climb  Search 

"Real  Time  A*  Search 


/*   Node  record  for  terrain  */ 


struct  node_rec 
{ 

int  xi; 

int  yi; 

float  x; 

float  y; 

int  state 


/*  for  graphics  */ 
/*  for  graphics  */ 

/*  UNTOUCHED,  OBSTACLE,  VISITED,  or  FRONTIER 


/* 
/* 


int  back_track_state; 

float  subtotal; 

float  dist_f rom_start ; 

float  dist_f rom_goal; 

float  dist_f rom_current ; 

struct  node_rec  *predecessor; 

int  f rontier_index; 

struct  node_rec  *qnext; 

struct  node_rec  *qreset; 


*/ 
*/ 


UNTOUCHED,  VISITED,  or  OBSTACLE 
includes  all  but  dist  from  current  &  move  away  */ 


/*  points  to  predecessor  for  a-star  search  */ 
/*  position  in  frontier  heap  */ 
/*  next  in  q  for  dist  from  current  DFS  */ 
/*  reset  link  list  after  back  DFS  */ 
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/*   Individual  record  for  population  */ 


}; 


struct  f actor_struct  { 

unsigned  int  place_holder  :  4 

unsigned  int  start_dist  :  4 

unsigned  int  goal_dist  :  4 

unsigned  int  current_dist  :  4 

unsigned  int  crowd_sides  :  4 

unsigned  int  crowd_diag  :  4 

unsigned  int  move_away  :  4 

unsigned  int  momentum  :  4 

}  ; 

union  chrom_union  { 

struct  f actor_struct  factor; 
unsigned  int  alleles; 

}  ; 

struct  individual_struct  { 

union  chrom_union  chrom; 

int  fitness; 

float  fit_sum; 

int  previous_index; 
}  ; 


/*  global  variables  */ 
extern  int  heap_size; 

/*   Functions  listed  under  file  */ 

/*  astar.c  */ 

float  a_star ()  ; 

int  update_astar_f rontier ()  ; 

/*  rt_astar.c  */ 

float  rta_star(); 

int  update_rtastar_adjacent () ; 

int  insert  ()  ; 

/*  hill_climb.c  */ 

float  hill_climb() ; 

struct  node_rec  *move_adjacent 0 ; 

struct  node_rec  *f ind_best ( ) ; 

/*  test.c  */ 
int  test () ; 
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/*  train. c  */ 
int  train  ()  ; 

/*  tsearch.c  */ 
float  search () ; 

/*  psearch.c  */ 

float  psearch () ; 

int  p_update_f rontier_list ( )  ; 

struct  node_rec  *p_pick_best_f rontier () ; 

/*  bfsearch.c  */ 

float  bf search (); 

int  bf_update_f rontier_list () ; 

struct  node_rec  *bf _j3ick_best_f rontier () ; 

/*  tsetup.c  */ 

int  get_seed( ) ; 

unsigned  int  get_pers_chrom( ) ; 

int  read_density_f ile ( ) ; 

int  make_array () ; 

int  make_node ( ) ; 

struct  node_rec  *f ind_node () ; 

/*  tpopulation  */ 

int  create_population ( ) ; 

struct  individual_struct  *new_individual ( ) ; 

int  get_population () ; 

/*  tprint.c  */ 

int  print_density () ; 

int  print_node ( ) ; 

int  print_population ()  ; 

/*  tfrontier.c  */ 

int  update_f rontier_list ()  ; 

struct  node_rec  *pick_best_f rontier () ; 

int  update_adjacent_obstacles  0 ; 

float  update_dist () ; 

int  update_crowd_sides () ; 

int  update_crowd_diag ( ) ; 

int  calc_move_away ( ) ; 

int  calc_momentum() ; 

float  compute_subtotal ( ) ; 
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/*  theap.c  */ 
int  insert_heap ( ) ; 
int  delete_heap () ; 
int  move_heap ( ) ; 
int  swap ( ) ; 

/*  tmisc.c  */ 

float  compute_shortest ()  ; 

int  equalf ( )  ; 

/*  update_dist_start ()  ;  */ 

int  gen_xi () ; 

int  gen_yi ( ) ; 

int  show_least_nodes ( ) ; 

/*  evolve. c  */ 

int  evolve () ; 

int  create_mate_heap () ; 

int  crossover  (); 

int  allele_crossover () ; 

int  bit_crossover () ; 

int  get_mask(); 

/*  eheap.c  */ 
int  insert_mate_heap () / 
int  pop_mate_heap ( ) ; 
int  move_mate_heap ( ) ; 
int  swap_num ( ) ; 

/*tdisplay.c  */ 
int  initialize () ; 
int  draw_terrain ( ) ; 
int  show_mouse ( ) ; 
int  draw_nodes () ; 
int  draw_grid(); 
void  squaref () ; 
void  square  ()  ; 
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main.c 

/* 

File: 

main.c 

Programmer: 

g.b.  parker 

Environment : 

any 

Language : 

C 

Date: 

9  july  92 

Revised: 

Comments : 

This  is  the 

This  is  the  control  for  user  input  and  call  of  train  or  test. 
Ten  command  live  arguments  are  optional.   The  syntax  for  a  call  is  as  follows 
(the  0  argument  is  the  program  call) :   (0)  t;   (1)  0  for  train,  1  for  test; 
(2)  random_seed;  (3)  input  population  file  name;  (4)  input  terrain  density  file 
name;  (5)  start  region  on  X  axis;  (6)  start  region  on  Y  axis;  (7)  goal  region 
on  X  axis;  (8)  goal  region  on  Y  axis;  (9)  number  of  generations  if  training, 
iterations  if  testing;  (10)  cycles  per  generation  if  training,  array  position 
of  best  individual  in  the  GA_produced  population  for  testing;  (11)  file  name 
for  out  population  if  training,  hexadecimal  representation  of  best  Persistence 
search  scheme. 
*/ 

♦include  "ga_search.h" 


/*  *****************  rnain  *******************************************  */ 


main(  argc,  argv  ) 
int  argc; 
char  *argv [] ; 

{ 

struct  individual_struct  *individual [32]  ; 

int  arg_seed  =  0; 

char  arg_population [32] ; 

char  arg_population_out [32] ; 

char  arg_density [32]  ; 

int  choice  =  0; 

int  sx  =  2, 

int  sy  =  2, 

int  gx  =  4, 

int  gy  =  2, 

int  iterations  =  25; 

int  generations  =  3; 

int  cycles_per_generation  =  2; 

int  best_indi vidua 1  =  -1; 

unsigned  int  pers_chrom; 

char  hname [64] ; 

int  hnlength; 
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gethostname (  hname,  hnlength  ) ; 

strcpy(  arg_population, 

strcpy(  arg_population_out,  "popx.out 

strcpy(  arg_density,         " 


"  ) 
"  ) 
"  ) 


switch  (  argc  ) 
{ 

case  12: 

sscanf 

sscanf 
case  11: 

sscanf 

sscanf 
case  10: 

sscanf 

sscanf 
case  9: 

sscanf 
case  8 : 

sscanf 
case  7: 

sscanf 
case  6: 

sscanf 
case  5 : 

sscanf 
case  4: 

sscanf 
case  3: 

sscanf 


argvfll],  "%s",  arg_jpopulation_out  ); 
argv[ll],  "%x" ,  &pers_chrom  ); 

argv[10],  "%d",  &cycles_per_generation  ); 
argv[10],  "%d",  &best_individual  ); 

argv[9],  "%d",  Sgenerations  )/ 
argv[9],  "%d",  Siterations  ); 


argv[8],  "%d",  &gy  ) 

argv[7],  "%d",  &gx  ) 

argv[6],  "%d",  &sy  ) 

argv[5],  "%d",  &sx  ) 

argv[4],  "%s",  arg_density  ); 

argv[3],  "%s",  arg_population  ); 

argv[2],  "%d",  &arg_seed  ); 
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case  2 : 

33canf(  argv[l],  "%d",  fichoice  ); 
switch  (  choice  ) 
{ 
case  0 : 

train (individual, arg_seed, arg_population,  arg_density, sx, sy, gx,gy, 

generations, cycles_per_generation, arg_jpopulation_out, hname    ) ; 
put_population  (    individual,    arg_population_out    ); 
break; 
case   1 : 

individual [0]  =  NULL; 

test (individual, arg_seed, arg_population, arg_density, sx, sy, gx, gy, 

iterations, best_individual,pers_chrom) ; 
break; 
case  2 : 

train (individual, arg_seed, arg_jpopulation, arg_density, sx, sy, gx, gy, 

generations,  cycles_per_generation, arg_j?opulation_out, hname) ; 
test (individual, arg_seed, arg_population,  arg_density,  sx,  sy,  gx, gy, 

iterations, best_individual, pers_chrom) ; 
put_population (  individual,  arg_population_out  ); 
break; 
} 
break; 
case  1: 

train (individual, arg_seed, arg_population, arg_density, sx, sy, gx, gy, 

generations, cycles_per_generation, arg_j30pulation_out, hname) ; 
put_population (  individual,  arg_population_out  ); 
} 


57 


train.c 

/* 

File:  train.c 

Programmer:  g.b.  parker 

Environment :  any 

Language :  C 

Date:         9  July  92 

Revised: 

Comments:  Called  by  main  to  train  a  population  of  32  individuals.   If  no 
input  population,  a  random  one  is  generated. 
*/ 

♦include  "ga_search.h" 


/*  *****************  train  ******************************************  */ 

train (  individual,  arg_seed,  arg_population,  arg_density,  sx,  sy,  gx,  gy,  gen- 
erations, cycles_per_generation,  arg_population_out,  hname  ) 
struct  individual_struct  *individual [32] ; 
int  arg_seed; 
char  arg_population [32] ; 
char  arg_density [32]  ; 
int  sx,  sy,  gx,  gy; 

int  generations,  cycles_per_generation; 
char  arg_population_out [32] ; 
char  hname [64] ; 
{ 

int  density [16]  [16] ; 

struct  node_rec  *node [ 66] [ 66] ; 

int  gen,  cycle,  i,  rs,  short_count; 

float  shortest_path; 

int  dummy; 

rs  =  arg_seed  ?  arg_seed  :  get_seed ( ) ; 
srandom(rs);   /*  seed  the  random  generator  */ 

print f ("\nRandom  seed  is  %d",  rs  ); 

get_population (  individual,  arg_population  ); 

read_density_f ile (  density,  arg_density  ); 
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for(  gen  -  1;  gen  <=  generations;  gen++  )  { 
printf("\n  gen  =  %d    (cycle, rs)  =  ",  gen); 

for(  cycle  =  1;  cycle  <=  cycles_per_generation;  cycle++  )  { 
rs  =  rs  +  1 ; 
short_count  =  0; 

while (  (shortest_path  =  a_star (sx, sy, gx, gy,  rs,  density, node  )) 
>  (  BIG_NUMBER  -  1.0  )  )  { 
if(  short_count  >  1000000000  )  { 

print f ("\nPROGRAM  ABORTED  -  iteration  %d  -  no  shortest  path\n",i) 
return  (F) ; 
} 
else 

rs  =  rs  +  1; 
} 

printfC  (%d,  %d)",  cycle,  rs  ); 
for(  i=0;  i<32;  i++  )  { 
if (  cycle  ==  1  ) 

individual [i] ->fit_sum  =  shortest_path  /  search (  sx,sy,gx,gy, 
individual [i] ->chrom. factor, rs, density, node, Sdummy  ) ; 
else 

individual [i] ->fit_sum  =  individual [i] ->fit_sum  +  shortest_path  / 
search  (  sx, sy, gx, gy, individual [i] ->chrom. factor, rs, density, 
node, & dummy  ) ; 
if (  cycle  ==  cycles_per_generation  ) 

individual [i] ->fitness  =  (int) ( (individual [i] ->fit_sum  / 
cycles_per_generation)  *  100.0); 
} 
} 
if  (gen  ==  generations  -  1  ) 

cycles_per_generation  =  cycles_per_generation  +  10; 

evolve (  individual,  rs  ) ; 
put_rs (  rs  ) ; 

if  (  (  (gen  %  10)  ==  0  )  II  (gen  <  10)  ) 
put_gen (gen, arg_population_out, hname, rs) ; 

/*  put  gen  to  a  standard  update  file  */ 
if  (  (gen  %  50)  ==  0  )  { 

put_population  (    individual,    arg_j5opulation_out    ); 
} 
} 
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/*  *****************  put  gen  ****************************************  */ 
/*  Called  by  train  to  continually  store  status  information  to  a  file  in  the 
directory  of  execution  */ 

put_gen  (  gen,  arg_j?opulation_out,  hname,  rs  ) 

int  gen; 

char   arg_population_out [32] ; 

char   hname [64] ; 

int    rs ; 

{ 

FILE    *gen_file,    *fopen(); 

gen_file  =  f open ("running. update",  "a"); 

fprintf (gen_f ile, "  %s  %s  gen  =  %d   rs  =  %d\n",  hname,  arg_population_out, 
gen,  rs) ; 

f close (gen_f ile) ; 

} 


/*  *****************  put  rs  *****************************************  */ 
/*  Puts  random  seed  info  to  a  file  in  the  directory  of  execution  */ 

put_rs (  rs  ) 

int  rs; 

{ 

FILE  *rs_file,  *fopen(); 

rs_file  =  f open ("rs .update",  "a"); 

fprintf (rs_f ile, "  %d  ",  rs)  ; 

fclose (rs  file) ; 
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teste 

/* 

File:  test.c 

Programmer:   g.b.  parker 

Environment :  any 

Language :     C 

Date:  9  July  92 

Revised: 

Comments:     Called  by  main  to  perform  a  comparative  test  of  search 
strategies.   The  default  is  for  all  individuals  of  the  population  to  be 
tested,  unless  a  specific  individual  is  specified. 
*/ 

♦include  "ga_search .h" 
♦include  <sys/time.h> 


/*  *****************  test  *******************************************  */ 

test (  individual,  arg_seed,  arg_population,  arg_density,  sx,  sy,  gx,  gy, 

iterations,  best_individual,  pers_chrom  ) 
struct  individual_struct  *individual [32] ; 
int  arg_seed; 
char  arg_population [32] ; 
char  arg_density [32] ; 
int  sx,  sy,  gx,  gy; 
int  iterations; 
int  best_individual; 
unsigned  int  pers_chrom; 
I 

struct  individual_struct  *best_first; 

struct  individual_struct  *persistence_search; 

int  density [16] [16] ; 

struct  node_rec  *node [ 66]  [66]  ; 

int  i,  rs,  k,  short_count; 

float  shortest_path; 
float  realtime_astar_f it_sum; 
int  realtime_astar_f itness; 
float  hill_climb_f it_sum; 
int  hill  climb  fitness; 
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float  temp; 
float  ga_t  =  0.0; 
float  ga_ticks  =  0.0; 
float  bf_t  =  0.0; 
float  bf_ticks  =  0.0; 
float  pers_t  =  0.0; 
float  pers_ticks  =  0.0; 
float  hc_t  =  0.0; 
float  hc_ticks  =  0.0; 
float  rta_t  =  0.0; 
float  rta_ticks=  0.0; 

long  sec,  usee; 

struct  timeval  *tvp  ■  (struct  timeval  *) malloc (sizeof (struct  timeval) ) ; 
struct  timezone  *tzp  =  (struct  timezone  *) malloc (sizeof (struct  timezone) ) ; 

best_first  = 

(struct  individual_struct  *)  malloc (sizeof (struct  individual_struct) ) ;; 
persistence_search  = 

(struct  individual_struct  *) malloc (sizeof (struct  individual_struct) ) ; 
best_first->chrom. alleles     =  OxeOlOOOOO; 
persistence_search->chrom. alleles  =  pers_chrom; 
best_f irst->f it_sum     =  0.0; 
persistence_search->f it_sum  =  0.0; 

rs  =  arg_seed  ?  arg_seed  :  get_seed(); 
srandom(rs);   /*  seed  the  random  generator  */ 

printf ("\nRandom  seed  is  %d",  rs)  ; 
if  (  individual [0]  ==  NULL  ) 

get_population (  individual,  arg_population  ); 

for(k=0;k<32;k++) 

individual [k] ->fit_sum  =  0.0; 

read_density_f ile (  density,  arg_density  ) ; 

printf ("\n (iteration,  rs)  ")  ; 
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for(  i  =  1;  i  <=  iterations;  i++  )  ( 
rs  =  rs  +  1; 
short_count  =  0; 

while  (  (shortest_path  =  a_star(  sx,  sy,  gx,  gy,  rs,  density,  node  )) 
>  (  BIG_NUMBER  -  1.0  )  )  { 
if(  short_count  >  1000000000  )  { 

printf("\n  PROGRAM  ABORTED  -  iteration  %d  -  no  shortest  path\n",  i  ); 
return (F) ; 
} 
else  { 

printf ("  (%d,  %d)  ",i,rs  ); 
rs  =  rs  +  1; 
} 
} 

printf ("  (%d,  %d)  ",i,rs  ); 
if  (best_individual  ==  -1) 
for(k=0;k<32;k++)  { 

temp  =  search (  sx, sy, gx, gy, individual [k] ->chrom. factor, rs, density, 

node, &ga_t  ) ; 
individual [k] ->fit_sum  = 

individual [k] ->fit_sum  +  shortest_path  /  temp; 
if  (temp  <  shortest_path  ) 

printf ("\nSHORTEST  PATH  >  ACTUAL  PATH  " ) ; 
I 
else  { 

k  =  best_individual; 

temp  =  search (  sx, sy, gx, gy, individual [k] ->chrom. factor, rs, density, 

node, &ga_t  ) ; 
individual [k] ->fit_sum  =  individual [k] ->fit_sum  +  shortest_path  /  temp; 
if  (temp  <  shortest_path  ) 

printf ("\nSHORTEST  PATH  >  ACTUAL  PATH  ") ; 
ga_ticks  =  ga_ticks  +  ga_t; 
best_f irst->f it_sum  =  best_f irst->f it_sum  +  shortest_path  / 

bf search  (  sx, sy, gx, gy, best_f irst->chrom. factor, 
rs, density, node, &bf_t  ); 
bf_ticks  =  bf_ticks  +  bf_t; 

persistence_search->f it_sum  =  persistence_search->f it_sum  + 
shortest_path  /  psearch(  sx,sy,gx,gy, 

persistence_search->chrom. factor,    rs, density, node, &pers_t    ); 
pers_ticks   =  pers_ticks   +  pers_t; 
hill_climb_f it_sum  =   hill_climb_f it_sum  +   shortest_jDath   / 

hill_climb(  sx,  sy,  gx,  gy,  rs,  density,  node,  &hc_t  ); 
hc_ticks  =  hc_ticks  +  hc_t; 
realtime_astar_f it_sum  =  realtime_astar_f it_sum  +  shortest_path  / 

rta_star (  sx,  sy,  gx,  gy,  rs,  density,  node,  &rta_t  ); 
rta_ticks  =  rta_ticks  +  rta_t; 
} 
} 
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if  (best_individual  ==  -1)  /*  no  best  individual  input  */ 
for(k=0;k<32;k++) 

individual [k] ->fitness  = 

(int) ( (individual [k]->fit_sum  /  iterations)  *  100.0); 
else 

individual [k] ->fitness  = 

(int) ( (individual [k] ->fit_sum  /  iterations)  *  100.0); 

best_first->fitness  =  (int) ( (best_first->f it_sum  /  iterations)  *  100.0); 
persistence_search->f itness  = 

(int) ( (persistence_search->fit_sum  /  iterations)  *  100.0); 
realtime_astar_f itness  = 

(int) ( (realtime_astar_f it_sum  /  iterations)  *  100.0); 
hill_climb_f itness  =  (int) ( (hill_climb_f it_sum  /  iterations)  *  100.0); 

if  (best_individual  ==  -1) 

print_population (  individual  ); 
else  ( 

printf("\n  ga-produced  %2d  %7.3f  %7.3f    %x",  individual [k] ->f itness, 

individual [k] ->f it_sum,  ga_ticks,  individual [k] ->chrom. alleles) ; 
printf("\n  best_first   %2d  %7.3f  %7.3f    %x",  best_first->f itness, 

best_f irst->f it_sum,  bf_ticks,  best_first->chrom. alleles) ; 
printf ("\n  persistence  %2d  %7.3f  %7.3f    %x", 

pers is tence_search->f itness,  persistence_search->f it_sum, 
pers_ticks,  persistence_search->chrom. alleles) ; 
printf ("\n  hill  climb   %2d  %7.3f  %7.3f",  hill_climb_f itness, 

hill_climb_f it_sum,  hc_ticks  ) ; 
printf ("\n  RTA  star    %2d  %7.3f  %7.3f",  realtime_astar_f itness, 
realtime_astar_f it_sum,  rta_ticks  ) ; 
} 
printf ("\n") ; 
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***********  * 


tsetup.c 

/* 

File:  tsetup.c 

Programmer:  g.b.  parker 

Environment :  any 

Language :  C 
Date:         9  July  92 
Revised: 

Comments:  Setup  functions 
*/ 

♦include  "ga_search.h" 

/*  *****************  get  seed  **************************** 
/*  Interfaces  with  user  to  get  random  seed  */ 

get_seed() 
{ 

char  nl[l];    /*  absorbs  new_line  after  seed  entry  */ 

int  rand_seed; 

printf ("\nEnter  random  seed  or  0  (system  assign  seed):  ") ; 

scanf ("%d", &rand_seed) ; 

gets (nl) ; 

if  (rand_seed  !=  0) 

return  rand_seed; 
else 

return  getpidO; 


/*  *****************  qet  pers  chrom  *********************************  */ 
/*  Interfaces  with  user  to  get  Persistence  chromosome  */ 

unsigned  int  get_pers_chrom( ) 
( 

char  nl[l];    /*  absorbs  new_line  after  seed  entry  */ 

unsigned  int  pers_chrom; 

printf  ("\nEnter  Persistence  chromosome  (8  hex  digits)  or  0  (eOHOOOO):  ")  ; 

scanf  ("%x", &pers_chrom) ; 

gets (nl) ; 

if  (pers_chrom  !=  0) 

return (  pers_chrom  ) ; 
else 

return  (  OxeOllOOOO  ); 
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/*  *****************  read  density  file  ******************************  */ 
/*  Reads  density  file  from  execution  directory  */ 

read_density_f ile (  density,  file_name  ) 
int  density[16] [16] ; 
char  f ile_name [32] ; 

( 

int  i,  j; 

FILE  *density_f ile,  *fopen(); 

int  not_end  =  T; 

int  node_density  =  STANDARD_DENSITY; 

if(  file_name[0]  ==  *  ■)  { 

printf ("\nEnter  density  file  name:  "); 

gets (file_name) ; 
} 

if  (  (density_f ile  =  fopen (f ile_name, "r" ) )  ==  NULL  )  { 

printf ("\nThe  file  does  not  exist,  standard  densities  being  used.")/ 
not_end  =  F; 
} 

for  (  j  =  15  ;  j>=0  ;  j—  )  { 
for  (  i=0  ;  i<=15  ;  i++  )  { 

if  (  not_end  &&  (  f scanf (density_f ile, "%x" ,  &node_density)  !=  EOF  )  ) 

density [i] [ j]  =  node_density; 
else  { 

density [i] [j]  =  STANDARD_DENSITY; 
not_end  =  F; 
} 
} 
} 

if (density_file  !=  NULL) 
fclose (density_f ile) ; 
} 
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/*  *****************  make  array  *************************************  */ 
/*  Creates  the  node  array  on  initial  use,  then  resets  records  after  that  */ 

make_array(  density,  node  ) 

int  density [16] [16] ; 

struct  node_rec  *node [ 66]  [  66] ; 

I 

static  int  first  =  T;   /*  indicates  if  first  time  to  make  array  */ 

int  i,  j; 

for(i=0;i<=65;i++)  { 

make_node (  node,  i,  0,  first  ); 
make_node (  node,  i,  65,  first  ) ; 
node[i] [0]->state  =  OBSTACLE; 
nodefi] [65]->state  =  OBSTACLE; 
} 

for( j=l; j<=64; j++)  { 

make_node (  node,  0,  j,  first  ); 
make_node (  node,  65,  j,  first  ); 
node[0] [ j] ->state  =  OBSTACLE; 
node[65] [ j] ->state  =  OBSTACLE; 
} 

for(j=l; j<=64; j++) 

for(i=l;i<=64;i++)  { 

make_node (  node,  i,  j,  first  ); 
if  (randl6()  <  density [ (i-1) /4]  [(  j-1) /4] ) 
node [i] [j]->st ate  =  OBSTACLE; 
} 
first  =  F; 
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/*  *****************  make  node  **************************************  */ 
/*  Resets  single  node  information  */ 

make_node (  node,  xi,  yi,  first  ) 

struct  node_rec  *node [66] [ 66] ; 

int  xi ; 

int  yi; 

int  first;   /*  T  of  F  */ 

{ 

int  k; 

if(  first  )   { 

node[xi][yi]  =  (struct  node_rec  *) malloc (sizeof (struct  node_rec) ) ; 

node[xi] [yi]->xi  =  xi; 

node [xi] [yi] ->yi  =  yi; 

node [xi] [yi] ->x  =  (float)xi; 

node [xi] [yi] ->y  =  (float)yi; 


} 


/*  OBSTACLE,  VISITED,  or  FRONTIER  */ 
/*  VISITED,  or  FRONTIER  */ 


nodefxi] [yi]->state  =  UNTOUCHED; 

node [xi] [yi] ->back_track_state  =  UNTOUCHED; 

node [xi] [yi] ->subtotal  =  BIG_NUMBER; 

node [xi] [yi] ->dist_f rom_start  -  BIG_NUMBER; 

node [xi] [yi] ->dist_f rom_goal  =  BIG_NUMBER; 

node [xi] [yi] ->dist_f rom_current  =  0.0; 

node [xi] [yi] ->predecessor  =  node [xi] [yi] ; 

/*  points  to  predecessor  for  a-star  search  */ 
nodefxi] [yi] ->f rontier_index  =  NA;    /*  not  have  index  to  f rontier_heap  */ 
node [xi] [yi] ->qnext  =  NULL; 

/*  points  to  next  in  q  for  dist  from  current  DFS  */ 
node [xi] [yi] ->qreset  =  NULL;    /*  reset  link  list  after  back  DFS  */ 
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/*  *****************  find  node  **************************************  */ 
/*  Picks  a  random  node  in  the  designated  density  area.  Used  to  identify 
start  and  goal  nodes.  */ 

struct  node_rec  *find_node(  node,  dens_col,  dens_row) 
struct  node_rec  *node [66]  [66]  ; 
int  dens_col; 
int  dens_row; 

{ 

int  k,  xi,  yi,  base_x,  base_y; 

base_x  =  (dens_col  *  4)  +  1; 
base_y  =  (dens_row  *  4)  +  1; 

for  (k=0;k<100;k++)  { 

xi  =  base_x  +  (randl6()  %  4) ; 
yi  =  basely  +  (randl6()  %  4) ; 

if  (node[xi] [yi]->state  ==  UNTOUCHED) 
return  node[xi] [yi] ; 

} 

return  node[xi] [yi] ; 
} 
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t  populations 

/* 

File:  tpopulation 

Programmer:  g.b.  parker 

Environment :  any 

Language :  C 

Date:  9  July  92 

Revised: 

Comments:  Functions  dealing  with  population  creation/storage 
*/ 

#include  "ga_search .h" 


/*  *****************  create  population  ******************************  */ 
/*  Generates  a  population  of  random  individuals  */ 

create_population  (  individual  ) 

struct  individual_struct  *individual [32] ; 

{ 

int  k; 

for  (k=0;k<32;k++)  { 

individual [k]  = 

(struct  individual_struct  *)malloc (sizeof (struct  individual_struct) ) ; 

new_individual (  individual [k]  ); 
} 


/*  *****************  new  individual  *********************************  */ 
/*  Sets  initial  values  of  individual  records  fields  */ 

struct  individual_struct  *new_individual (  ind  ) 

struct  individual_struct  *ind; 

{ 

ind->chrom. factor .place_holder  =  Oxf; 

ind->chrom. factor . st art_dist    =  randl6() 

ind->chrom. factor .goal_dist     =  randl6() 

ind->chrom. factor .current_dist  =  randl6() 

ind->chrom. factor .crowd_sides   =  randl6() 

ind->chrom. factor . crowd_diag   =  randl6() 

ind->chrom. factor .move_away    =  randl6() 

ind->chrom. factor .momentum     =  randl6() 

ind->fitness      =  0; 

ind->fit_sum      =  0.0; 

ind->previous_index  =  99; 
) 
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/*  *****************  get  population  *********************************  */ 
/*  Reads  population  from  a  file  */ 

get_population (  individual,  file_name  ) 
struct  individual_struct  *individual [32] ; 
char  f ile_name [32] ; 
{ 

int  i ; 

FILE  *population_f ile,  *fopen(); 

int  not_end  =  T; 

unsigned  int  alleles; 

if    (    file_name[0]    —    •    •    )     { 

printf ("\nEnter  population  file  name:  ") ; 

gets (f ile_name) ; 
} 

if  (  (population_f ile  =  f open (f ile_name, "r") )  ==  NULL  )  { 

printf ("\nThe  file  does  not  exist,  random  population  being  used."); 
create_population (  individual  ); 

} 

else  { 

for  (  i=0  ;  i<32  ;  i++  )  { 
individual [i]  = 

(struct  individual_struct  *) malloc (sizeof (struct  individual_struct) ) ; 
individual [i] ->fitness      =  0; 
individual [i] ->fit_sum      =  0.0; 
individual [i] ->previous_index  =  99; 

if  (  not_end  &&  (  f scanf (population_file, "%x", &alleles)  !=  EOF  )  ) 
individual [i] ->chrom. alleles  =  alleles; 

else  { 

new_individual (  individual [i]  ); 
not_end  =  F; 
} 
I 

f close (population_f ile) ; 
} 
I 
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/*  *****************  put  population  *********************************  */ 
/*  Puts  the  population  to  a  designated  file  */ 

put_population (  individual,  file_name  ) 
struct  individual_struct  *individual [32] ; 
char  f ile_name [32] ; 
{ 

int  i  ; 

FILE  *population_f ile,  *fopen(); 

if(  file_name[0]  ==''){ 

print f ("\nEnter  output  population  file  name:  "); 
gets (f ile_name) ; 

} 

population_f ile  =  f open (f ile_name, "w") ; 

for  (  i=0  ;  i<32  ;  i++  ) 

fprintf (population_f ile, "%x\n", individual [i] ->chrom. alleles) ; 

f close (population_f ile) ; 
} 
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astar.c 

/* 

File:  astar.c 

Programmer:  g.b.  parker 

Environment:  any 

Language :  C 
Date:         9  July  92 
Revised: 

Comments:  A-star  search  -  Finds  the  shortest  path 
*/ 

♦include  "ga_search.h" 

int  heap  size; 


/*  *****************  a  star  ***************************************** 

float  a_star(  sx,  sy,  gx,  gy,  random_seed,  density,  node  ) 

int  sx,sy,gx,gy;     /*  position  in  density  array  for  start  &  goal  */ 

int  random_seed; 

int  density[16] [16] ; 

struct  node_rec  *node [ 66] [ 66] ; 

I 

struct  node_rec  *current,  *start,  *goal; 

struct  node_rec  *f rontier_heap [40  96] ; 

int  k; 

#ifdef  IRIS 

union  chrom_union  dummy_cu; 

dummy_cu .alleles  =  0; 
#endif 

heap_size  =  0; 

s random (random_seed) ;   /*  seed  the  random  generator  */ 

make_array (density,  node) ; 

start  =  f ind_node (node,  sx,  sy) ; 

start->state  =  START; 

start->dist_f rom_start  =  0.0; 

goal  =  f ind_node (node,  gx,  gy) ; 

goal->state  =  GOAL; 

for  (k=0;k<4  0  96;k++) 

f rontier_heap [k]  =  NULL; 
current  =  start; 
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#ifdef  IRIS 

/*  initialize  the  IRIS  system  */ 

initialize  (ASTITLE)  ; 

draw_terrain (node, start, goal, current , goal->dist_f rom_start, 
dummy_cu . factor) ; 
#endif 

while  (  current  !=  goal)  { 

update_astar_f rontier (  node, current , frontier_heap, start, goal  ) ; 

if (  heap_size  ==  0)  break; 

current  =  f rontier_heap [0] ; 

delete_heap(  f rontier_heap,  f rontier_heap [0]  ) ; 
} 

#ifdef  IRIS 

draw_terrain (node, start, goal, current, goal->dist_f rom_start, 
dummy_cu . factor) ; 

sleep (  SLEEPTIME  ) ; 
#endif 

goal->state  =  GOAL; 
start->state  =  START; 

return (  goal->dist_f rom_start  ); 
» 
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/*  *****************  update  astar  frontier  **************************  */ 

update_astar_frontier (  node,  c,  front ier_heap,  start,  goal  ) 

struct  node_rec  *node [ 66]  [ 66]  ; 

struct  node_rec  *c;         /*  current  */ 

struct  node_rec  *f rontier_heap [ 4096] ; 

struct  node_rec  *start; 

struct  node_rec  *goal; 

{ 

int  xi,  yi,  base_xi,  base_yi,  top_xi,  top_yi; 

base_xi  =  c->xi  ==  1  ?  1  :  c->xi  -  1 

base_yi  =  c->yi  ==  1  ?  1  :  c->yi  -  1 

top_xi  =  c->xi  ==  64  ?  64  :  c->xi  +  1 

top_yi  =  c->yi  ==  64  ?  64  :  c->yi  +  1 

for  (xi=base_xi;xi<=top_xi;xi++) 
for  (yi=base_yi;yi<=top_yi;yi++) 

if  (  (node [xi] [yi]->state  ==  UNTOUCHED)  || 
(node[xi] [yi]->state  ==  GOAL)  )  { 
node[xi] [yi] ->state  =  FRONTIER; 
node [xi] [yi] ->predecessor  =  c; 

node[xi] [yi] ->dist_f rom_goal   =  update_dist (node [xi] [yi],goal); 
node [xi] [yi] ->dist_f rom_start  = 

c->dist_f rom_start  +  update_dist (node [xi] [yi] , c) ; 
node [xi] [yi] ->subtotal  =  node [xi] [yi] ->dist_f rom_goal  + 

node[xi] [yi] ->dist_f rom_start ; 
insert_heap(  f rontier_heap,  node[xi][yi]  ); 
1 

else  if  (  node [xi] [yi] ->state  ==  FRONTIER  )  { 
if  (  node  [xi]  [yi]  ->dist_f  rorn_start  > 

(c->dist_f rom_start  +  update_dist (node [xi] [yi] , c) )  )  { 
node [xi] [yi] ->dist_f rom_start  = 

c->dist_f rom_start  +  update_dist (node [xi] [yi] , c) ; 
node [xi] [yi] ->subtotal  =  node [xi] [yi] ->dist_f rom_goal  + 

node[xi] [yi] ->dist_f rom_start ; 
move_heap (  f rontier_heap,  node [xi] [yi] ->f rontier_index  ); 
} 
} 
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hillclimbx 

/* 

File: 

Programmer: 

Environment 

Language : 

Date: 

Revised: 

Comments : 

tracks 

*/ 


hill_climb.c 
g.b.  parker 
any 
C 
9  July  92 

Hill  climb  search  -  Finds  best  adjacent  frontier  node  or  back- 


#include  "ga_search .h" 
♦include  <sys/time.h> 


/*  *****************  hill  climb  *************************************  */ 

float  hill_climb(  sx,  sy,  gx,  gy,  random_seed,  density,  node,  ticks  ) 

int  sx,sy,gx,gy;     /*  position  in  density  array  for  start  &  goal  */ 

int  random_seed; 

int  density [16]  [16]  ; 

struct  node_rec  *node [ 66]  [  66] ; 

float  *ticks; 

{ 

struct  node_rec  *current,  *next,  *start,  *goal; 

int  k  =  1; 

long  sec,  usee; 
static  struct  timeval  *tvp; 
static  struct  timezone  *tzp; 
static  int  first  =  T; 

#ifdef  IRIS 

union  chrom_union  dummy_cu; 

dummy_cu. alleles  =  0; 
#endif 

if(  first  )  { 

tvp  =  (struct  timeval  *) malloc (sizeof (struct  timeval)); 

tzp  =  (struct  timezone  *) malloc (sizeof (struct  timezone)); 

first  =  F; 
} 

srandom(random_seed) ;   /*  seed  the  random  generator  */ 
make_array (density,  node); 
start  =  f ind_node (node,  sx,  sy)  ; 
start->state  =  START; 
start->dist  from  start  =  0.0; 
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goal  =  f ind_node (node,  gx,  gy) ; 
goal->state  =  GOAL; 
current  -  start; 
current->state  =  CURRENT; 

#ifdef  IRIS 

/*  initialize  the  IRIS  system  */ 

initialize (HCTITLE) ; 
#endif 

gettimeofday (tvp,  tzp); 
sec  =  tvp->tv_sec; 
usee  =  tvp->tv_usec; 

while  (  current  !=  goal)  { 

next  =  move_adjacent (  node,  current,  start,  goal  ); 

if  (  next  !=  NULL  ) 

next->predecessor  =  current; 
else  if  (  current->predecessor  !=  NULL  ) 

next  =  current->predecessor; 
else 

printf ("\nNO  SOLUTION  -  hill  climb  search"); 

next->dist_f rom_start  = 

current->dist_f rom_start  +  update_dist (  current,  next  )  ; 
current->state  =  VISITED; 
current  =  next; 
current->state  =  CURRENT; 
#ifdef  IRIS 

draw_terrain (node, start, goal, current, next->dist_f rom_start, 
dummy _cu. factor) ; 
#endif 
} 

gettimeofday (tvp,  tzp) ; 

*ticks  =  (float) (tvp->tv_sec  -  sec)  +  (tvp->tv_usec  -  usee) /1000000 . 0; 

#ifdef  IRIS 

sleep (  SLEEPTIME  ) ; 
#endif 

goal->state  =  GOAL; 
start->state  =  START; 

return (  goal->dist_f rom_start  ); 
} 
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1 

:  c->xi 

- 

1 

1 

:  c->yi 

- 

1 

64 

:  c->xi 

+ 

1 

64 

:  c->yi 

+ 

1 

/•  *****************  move  adjacent  **********************************  */ 

struct  node_rec  *move_adjacent (  node,  c,  start,  goal  ) 

struct  node_rec  *node [66]  [66]  / 

struct  node_rec  *c;         /*  current  */ 

struct  node_rec  *start; 

struct  node_rec  *goal; 

{ 

struct  node_rec  *best; 

int  xi,  yi,  base_xi,  base_yi,  top_xi,  top_yi; 

best  =  NULL; 

base_xi  =  c->xi  ==  1  ? 

base_yi  =  c->yi  ==  1  ? 

top_xi   =  c->xi  ==  64  ? 

top_yi   =  c->yi  ==  64  ? 

for  (xi=base_xi;xi<=top_xi;xi++) 
for  (yi=base_yi;yi<=top_yi;yi++) 
switch  (  node [xi] [yi] ->state  ) 

{ 

case  UNTOUCHED: 
case  GOAL: 

node [xi] [yi] ->state  =  FRONTIER; 

node [xi] [yi] ->di3t_f rom_goal   =  update_dist (  node [xi] [yi] ,  goal   ); 

best  =  f ind_best (  best,  node[xi][yi]  ); 

break; 
case  FRONTIER: 

best   =   f ind_best (   best,    node[xi][yi]    ); 

break; 
I 
return (  best  ) ; 
} 


/*  *****************  find  best  **************************************  */ 
/*  Assigns  the  input  node  to  best  if  appropriate  */ 

struct  node_rec  *find_best(  best,  n  ) 
struct  node_rec  *best; 
struct  node_rec  *n; 
{ 

if  (  best  ==  NULL  ) 
best  =  n; 

else  if  (  n->dist_f rom_goal  <  best->dist_f rom_goal  ) 
best  =  n; 

return (  best  ) ; 
} 
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rt_astar.c 

/* 

File:  rt_astar.c 

Programmer:  g.b.  parker 

Environment :  any 

Language :  C 

Date:  20  feb  92 
Revised:       2  apr  92 

Comments:  RTA-star  search  -  Finds  best  adjacent  node  visited  or  frontier 
*/ 

♦include  "ga_search .h" 
♦include  <sys/time.h> 

int  heap_size; 


/*  *****•**•***•**•*  rta  star 


•a*************************************  * 


float  rta_star(  sx,  sy,  gx,  gy,  random_seed,  density,  node,  tick3  ) 

int  sx,sy,gx,gy;     /*  position  in  density  array  for  start  &  goal  */ 

int  random_seed; 

int  density [16] [16] ; 

struct  node_rec  *node [ 66] [66] ; 

float  *ticks; 

{ 

struct  node_rec  *current,  *start,  *goal; 

struct  node_rec  *best_two [2] ; 

int  k  =  1 ; 

long  sec,  usee; 

static  struct  timeval  *tvp; 

static  struct  timezone  *tzp; 

static  int  first  =  T; 
#ifdef  IRIS 

union  chrom_union  dummy_cu; 

dummy_cu .alleles  =  0; 
#endif 

if(  first  )  { 

tvp  =  (struct  timeval  *)malloc (sizeof (struct  timeval)); 

tzp  =  (struct  timezone  *) malloc (sizeof (struct  timezone)); 

first  =  F; 
) 

heap_size  =  0; 

s random (random_seed) ;   /*  seed  the  random  generator  */ 

make_array (density,  node); 
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start  =  f ind_node (node,  sx,  sy)  ; 
3tart->state  =  START; 
start->dist_f rom_start  =  0.0; 
goal  =  f ind_node (node,  gx,  gy) ; 
goal->state  =  GOAL; 
current  =  start; 
current->state  =  CURRENT; 

#ifdef  IRIS 

/*  initialize  the  IRIS  system  */ 

initialize (RTASTITLE)  ; 
#endif 

gettimeofday (tvp,  tzp); 
sec  =  tvp->tv_sec; 
usee  =  tvp->tv_usec; 

while  (  current  !=  goal)  { 
best_two[0]  =  NULL; 
best_two[l]  =  NULL; 

update_rtastar_adjacent (  node,  current,  best_two,  start,  goal  ); 
if (  best_two[0]  ==  NULL  )  break; 
if(  best_two[l]  ==  NULL  ) 

current->dist_from_goal  =  BIG_NUMBER; 
else 

current->dist_f rom_goal  =  best_two [1] ->subtotal; 
current->state  =  VISITED; 
best_two [0] ->dist_f rom_start  = 

current->dist_f rom_start  +  update_dist (current, best_two [0] ) 
current  =  best_two[0]; 
current->state  =  CURRENT; 

#ifdef  IRIS 

draw_terrain (node, start, goal, current, best_two [0] ->dist_f rom_start, 
dummy_cu . factor) ; 
#endif 
} 

gettimeofday (tvp,  tzp); 

*ticks  =  (float) (tvp->tv_sec  -  sec)  +  (tvp->tv_usec  -  usee) /1000000 .0; 

#ifdef  IRIS 

sleep (  SLEEPTIME  ) ; 
#endif 

goal->state  =  GOAL; 

start->state  =  START; 

return (  goal->dist_f rom_start  ); 

} 
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/*  *****************  update  rtastar  adjacent  ************************  * 

update_rta3tar_adjacent (  node,  c,  best_two,  start,  goal  ) 

struct  node_rec  *node [66] [66] ; 

struct  node_rec  *c;         /*  current  */ 

struct  node_rec  *best_two [2] ; 

struct  node_rec  *start; 

struct  node_rec  *goal; 

{ 

int  xi,  yi,  base_xi,  base_yi,  top_xi,  top_yi; 

base_xi  =  c->xi  ==  1  ?  1  :  c->xi  -  1 

base_yi  =  c->yi  ==  1  ?  1  :  c->yi  -  1 

top_xi  =  c->xi  ==  64  ?  64  :  c->xi  +  1 

top_yi  =  c->yi  ==  64  ?  64  :  c->yi  +  1 

for  (xi=base_xi;xi<=top_xi;xi++) 
for  (yi=base_yi;yi<=top_yi;yi++) 
switch  (  node [xi] [yi] ->state  ) 
{ 

case  UNTOUCHED: 
case  GOAL: 

node [xi] [yi] ->state  =  FRONTIER; 

node [xi] [yi] ->predecessor  =  c; 

node [xi] [yi] ->dist_f rom_goal   =  update_dist (  node [xi] [yi] ,  goal   )/ 

node [xi] [yi] ->subtotal  = 

node [xi] [yi] ->dist_f rom_goal  +  update_dist (  node [xi] [yi] ,  c  ); 
insert (  best_two,  node[xi][yi]  ) ; 
break; 
case  FRONTIER: 
case  VISITED: 

node [xi] [yi] ->subtotal  = 

node [xi] [yi] ->dist_f rom_goal  +  update_dist (  node [xi] [yi] ,  c  ); 
insert  (  best_two,  node[xi][yi]  ); 
break; 
case  CURRENT: 

break; 
} 
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/*  *****************  insert  *****************************************  */ 
/*  Assigns  the  input  node  to  best  or  second_best  as  appropriate  */ 

insert  (  best_two,  n  ) 

struct  node_rec  *best_two [2] ; 

struct  node_rec  *n; 

{ 

if  (  best_two[0]  ==  NULL  ) 

best_two[0]  =  n; 
else  if  (  n->subtotal  <  best_two [0] ->subtotal  )  { 
best_two[l]  =  best_two[0]; 
best_two[0]  =  n; 
} 
else  if  (  best_two[l]  —  NULL  ) 

best_two[l]  =  n; 
else  if  (  n->subtotal  <  best_two [1] ->subtotal  ) 
best  two[l]  =  n; 
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bfsearch.c 

/* 

File:        bfsearch.c 

Programmer:   g.b.  parker 

Environment :  any 

Language :     C 

Date:         9  July  92 

Revised: 

Comments:     Best  First  Search  -  modified  standard  bfs  for  use  with  real- 
time search 
*/ 

♦include  "ga_search .h" 
♦include  <sys/time.h> 


/*  *****************  bfsearch  ***************************************  */ 

float  bfsearch (  sx,  sy,  gx,  gy,  ind_chrom_f actor,  random_seed,  density,  node, 

ticks  ) 

int  sx,sy,gx,gy;     /*  position  in  density  array  for  start  &  goal  */ 

struct  factor_struct  ind_chrom_factor; 

int  random_seed; 

int  density [16]  [16]  ; 

struct  node_rec  *node [ 66]  [66]  ; 

float  *ticks; 

{ 

struct  node_rec  *current,  *previous,  *start,  *goal; 

struct  node_rec  *f rontier_heap [ 4096] ; 

long  sec,  usee; 
static  struct  timeval  *tvp; 
static  struct  timezone  *tzp; 
static  int  first  =  T; 

float  dist_traveled  =  0.0; 
int  k; 

#ifdef  IRIS 

union  chrom_union  dummy_cu; 

dummy_cu .alleles  =  0; 
#endif 

if(  first  )  { 

tvp  =  (struct  timeval  *) malloc (sizeof (struct  timeval)); 
tzp  =  (struct  timezone  *) malloc (sizeof (struct  timezone)); 
first  =  F; 

} 
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srandom(random_seed) ;   /*   seed  the  random  generator  */ 

heap_size  =  0; 

make_array (density,  node); 

start  =  f ind_node (node,  sx,  sy) ; 
start->state  =  START; 
start ->dist_from_st art  =  0.0; 
goal  =  f ind_node (node,  gx,  gy) ; 
goal->state  =  GOAL; 

for  (k=0;k<4096;k++) 

f rontier_heap [k]  =  NULL; 

current  =  start; 

previous  =  start; 
#ifdef  IRIS 

/*  initialize  the  IRIS  system  */ 

initialize (BFTITLE) ; 
#endif 

gettimeofday (tvp,  tzp); 

sec  =  tvp->tv_sec; 

usee  =  tvp->tv_usec; 

while  (  current  !=  goal)  { 

bf_update_f rontier_list (node, current, previous, front ier_heap, start, goal) ; 
if(  heap_size  ==0)  { 

printf ("\nENDING  SEARCH  BEFORE  GOAL  -  no  more  frontiers- 
break; 
} 

previous  =  current; 

current  =  bf_pick_best_f rontier (node, current , front ier_heap, goal) ; 
dist_t raveled  =  dist_t raveled  +  current->dist_f rom_current ; 
#ifdef  IRIS 

draw_terrain (node, start, goal, current, dist_traveled, dummy_cu . factor) ; 
#endif 
> 

gettimeofday (tvp,  tzp); 

*ticks  =  (float)  (tvp->tv_sec  -  sec)  +  (tvp->ty_usec  -  usee) /1000000 . 0; 

#ifdef  IRIS 

sleep  (  SLEEPTIME  ) ; 
#endif 

goal->state  =  GOAL; 
start->state  =  START; 

return (  dist_traveled  ); 
1 
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/*  *****************  bf  update  frontier  list  ************************  */ 

bf_update_f rontier_list (  node,  c,  p,  f rontier_heap,  start,  goal  ) 

struct  node_rec  *node [66] [66] ; 

struct  node_rec  *c;         /*  current  */ 

struct  node_rec  *p;         /*  previous  */ 

struct  node_rec  *f rontier_heap [4096] ; 

struct  node_rec  * start; 

struct  node_rec  *goal; 

{ 

int  xi,  yi,  base_xi,  base_yi,  top_xi,  top_yi; 

float  old_subtotal; 

base_xi  =  c->xi  -  1 
base_yi  =  c->yi  -  1 
top_xi  =  c->xi  +  1 
top_yi   =  c->yi  +  1 

for  (xi=base_xi;xi<=top_xi;xi++) 
for  (yi=base_yi;yi<=top_yi;yi++) 

if  (  (nodefxi] [yi]->state  —  UNTOUCHED)  ||  (node [xi] [yi] ->state  ==  GOAL) 


)  ( 


node [xi] [yi] ->state  =  FRONTIER; 

node [xi] [yi] ->predecessor  =  c; 

node [xi] [yi] ->dist_f rom_goal   =  update_dist (  node [xi] [yi] ,  goal  ); 

node [xi] [yi] ->subtotal  =  node [xi] [yi] ->dist_f rom_goal; 

insert_heap (  f rontier_heap,  node[xi][yi]  ); 


/*  *****************  bf  pick  best  frontier  **************************  */ 
/*  Finds  best  frontier  node,  returns  it  */ 

struct  node_rec  *bf_pick_best_f rontier (  node,  current,  f rontier_heap,  goal  ) 

struct  node_rec  *node [66] [ 66] ; 

struct  node_rec  *current; 

struct  node_rec  *f rontier_heap [4096] ; 

struct  node_rec  *goal; 

i 

struct  node_rec  *best_ptr,  *q,  *qend,  *qreset; 

float  node_cost  =  BIG_NUMBER; 

float  norm  =  1.0;  /*  (current->dist_f rom_goal  /  16.0);    normalize  factor  */ 

int  xi,  yi,  k; 

float  steps; 

int  done  =  F; 
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} 


best_ptr  -  f rontier_heap [0] ; 

q  =  current; 

qend  =  current; 

qreset  =  current; 

current ->back_track_state  =  VISITED; 

current->dist_f rom_current  =  0.0; 

while  (  !done  &&  (q  !=  NULL)  )  { 

steps  =  q->dist_f rom_current  +  1.0; 
for(k=0;k<8;k++)  { 

if  (k==4) 

steps  =   q->dist_f rom_current  +  SQRT2; 

xi  =  gen_xi (  k,  q->xi  ) ; 

yi  =  gen_yi (  k,  q->yi  ) ; 

if(   (  (nodefxi] [yi]->state  ==  VISITED)  || 
(node[xi] [yi]->state  ==  FRONTIER)  || 
(node[xi] [yi] ->state  ==  START)  )  && 
node[xi] [yi] ->back_track_state  ==  UNTOUCHED  )  { 
node [xi] [yi] ->dist_f rom_current  =  steps; 
node [xi] [yi] ->back_track_state  =  VISITED; 
node [xi] [yi] ->qreset  =  qreset; 
qreset  =  node [xi] [yi] ; 

if  (  (node [xi] [yi]->state  ==  VISITED)  || 
(node [xi] [yi] ->state  ==  START)  )  { 
qend->qnext  =  node [xi] [yi] ; 
qend  =  node[xi] [yi] ; 
1 

else  (     /*  nodefxi] [yi]->state  ==  FRONTIER  */ 
if  (  node[xi] [yi]  ==  best_ptr  ) 
done  =  T; 
} 
} 
} 

q  =  q->qnext; 
} 

reset_back_track_state (  qreset  ) ; 
best_ptr->state  =  VISITED; 
delete_heap(  f rontier_heap,  best_ptr  ); 

return (  best_ptr  ); 
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Ipsearch.c 
/* 

File:  psearch.c 

Programmer:  g.b.  parker 

Environment:  any 

Language :  C 

Date:         9  July  92 

Revised: 

Comments:  Persistence  Search  -  uses  distance  to  goal  and  distance  to  cur- 
rent 

to  determine  best  frontier. 
*/ 

♦include  "ga_search .h" 
♦include  <sys/time.h> 


/*  *****************  psearch  ****************************************  */ 

float  psearch  (  sx,  sy,  gx,  gy,  ind_chrom_f actor,  random_seed,  density,  node, 

ticks  ) 

int  sx,sy,gx,gy;     /*  position  in  density  array  for  start  &  goal  */ 

struct  f actor_struct  ind_chrom_f actor; 

int  random_seed; 

int  density [16]  [16]  ; 

struct  node_rec  *node [66]  [ 66]  ; 

float  *ticks; 

( 

struct  node_rec  *current,  *previous,  *start,  *goal; 

struct  node_rec  *f rontier_heap [4096] ; 

long  sec,  usee; 
static  struct  timeval  *tvp; 
static  struct  timezone  *tzp; 
static  int  first  =  T; 
float  dist_traveled  =  0.0; 
int  k; 

if(  first  )  { 

tvp  =  (struct  timeval  *)malloc (sizeof (struct  timeval)); 
tzp  =  (struct  timezone  *)malloc (sizeof (struct  timezone)); 
first  =  F; 

} 

srandom(random_seed) ;   /*   seed  the  random  generator  */ 

heap_size  =  0; 

make_array (density,  node); 

start  =  f ind_node (node,  sx,  sy) ; 

start->state  =  START; 
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start ->dist_f rom_start  =  0.0; 
goal  =  f ind_node (node,  gx,  gy) ; 
goal->state  =  GOAL; 

for  (k=0;k<4096;k++) 

f rontier_heap [k]  =  NULL; 

current  =  start; 

previous  =  start; 
#ifdef  IRIS 

/*  initialize  the  IRIS  system  */ 

initialize (PERSTITLE) ; 
#endif 

gettimeofday (tvp,  tzp); 

sec  =  tvp->tv_sec; 

usee  =  tvp->tv_usec; 

while  (  current  !=  goal)  { 

if (  adjacent (current,  goal)  )  { 

dist_traveled  =  dist_traveled  +  update_dist (  current,  goal  ); 
current  =  goal; 
} 
else  ( 

p_update_f rontier_list (  node,  current,  previous,  f rontier_heap,  start, 

goal,  ind_chrom_factor  ) ; 
if (  heap_size  ==0)  { 

printf ("\nENDING  SEARCH  BEFORE  GOAL  %d  %d   -  no  more  frontier", 

current->xi,  current->yi) ; 
break; 
} 

previous  =  current; 
current  =  p_pick_best_f rontier (  node,  current,  f rontier_heap, 

goal,  ind_chrom_factor  ) ; 
dist_t raveled  =  dist_traveled  +  current->dist_f rom_current ; 
} 
#ifdef  IRIS 

draw_terrain (node, start, goal, current, dist_traveled, ind_chrom_f actor) ; 
#endif 
} 

gettimeofday (tvp,  tzp); 

*ticks  =  (float) (tvp->tv_sec  -  sec)  +  (tvp->tv_usec  -  usee) /1000000 . 0; 
#ifdef  IRIS 

sleep (  SLEEPTIME  ) ; 
#endif 

goal->state  =  GOAL; 

start->state  =  START; 

return (  dist_traveled  ); 

} 
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/*  *****************  p  update  frontier  list  *************************  */ 

p_update_frontier_list (  node,  c,  p,  f rontier_heap,  start,  goal,  factor  ) 

struct  node_rec  *node [ 66] [ 66] ; 

struct  node_rec  *c;         /*  current  */ 

struct  node_rec  *p;         /*  previous  */ 

struct  node_rec  *f rontier_heap [4096] ; 

struct  node_rec  *start; 

struct  node_rec  *goal; 

struct  f actor_struct  factor; 

{ 

int  xi,  yi,  base_xi,  base_yi,  top_xi,  top_yi; 

float  old_subtotal; 

base_xi  =  c->xi  -  1 

base_yi  =  c->yi  -  1 

top_xi   =  c->xi  +  1 

top_yi   =  c->yi  +  1 

for  (xi=base_xi;xi<=top_xi;xi++) 
for  (yi=base_yi;yi<=top_yi;yi++) 

if ( (node[xi] [yi]->state  ==  UNTOUCHED)  ||  (node [xi] [yi] ->state  ==  GOAL))  { 
nodefxi] [yi]->state  =  FRONTIER; 
node [xi] [yi] ->predecessor  -  c; 

node [xi] [yi] ->dist_f rom_goal   =  update_dist (  node [xi] [yi] ,  goal  ); 
node [xi] [yi] ->subtotal  = 

node [xi] [yi] ->dist_f rom_goal  *  f actor .goal_dist; 
insert_heap(  f rontier_heap,  nodefxi] [yi]  ); 
} 
} 

/*  *****************  p  pick  best  frontier  ***************************  */ 
/*  Finds  best  frontier  node,  returns  it  */ 

struct  node_rec  *p_jpick_best_f rontier (  node,  current,  f rontier_heap,  goal,  fac- 
tor ) 

struct  node_rec  *node [ 66] [66] ; 
struct  node_rec  *current; 
struct  node_rec  *f rontier_heap [4096] ; 
struct  node_rec  *goal; 
struct  factor_struct  factor; 
{ 

struct  node_rec  *best_ptr,  *q,  *qend,  *qreset; 

float  node_cost  =  BIG_NUMBER; 

float  norm  =  1.0;  /*  (current->dist_f rom_goal  /  16.0);    normalize  factor  */ 

float  lower_bound  =  f rontier_heap [0] ->subtotal; 

float  upper_bound  =  BIG_NUMBER; 

int  xi,  yi,  k; 

float  steps; 
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} 


best_ptr  -  current; 
q  =  current; 
qend  =  current; 
qreset  =  current; 

current->back_track_state  =  VISITED; 
current->dist_f rom_current  =  0.0; 

while  (  (lower_bound  <  upper_bound)  &&  (q  !=  NULL)  )  { 
steps  =  q->dist_f rom_current  +  1.0; 
for(k=0;k<8;k++)  { 
if  (k==4) 

steps  =   q->dist_f rom_current  +  SQRT2; 
xi  =  gen_xi (  k,  q->xi  ) ; 
yi  =  gen_yi (  k,  q->yi  ); 

if  (  (  (node[xi] [yi]->state  ==  VISITED)  || 
(node[xi] [yi]->state  ==  FRONTIER)  || 
(node [xi] [yi]->state  ==  START)  )  && 
nodefxi] [yi] ->back_track_state  ==  UNTOUCHED  )  { 
node [xi] [yi] ->dist_f rom_current  -   steps; 
node [xi] [yi] ->back_track_state  =  VISITED; 
node [xi] [yi] ->qreset  =  qreset; 
qreset  =  node [xi]  [yi]  ; 

if  (  (node[xi] [yi]->state  ==  VISITED)  || 
(nodefxi] [yi] ->state  ==  START)  )  { 
qend->qnext  =  node [xi] [yi] ; 
qend  =  node[xi] [yi]  ; 
} 

else  {     /*  node [xi] [yi] ->state  ==  FRONTIER  */ 
node_cost  =  node [xi] [yi] ->subtotal  + 

node [xi] [yi] ->dist_f rom_current  *  factor .current_dist; 

if  (node_cost  <  upper_bound)  { 
upper_bound  =  node_cost; 
best_ptr  =  node [xi]  [yi]  ; 
} 
} 
} 
} 

q  =  q->qnext; 

lower_bound  =  front ier_heap [0] ->subtotal  +  steps  *  factor .current_dist; 
} 

reset_back_track_state (  qreset  ); 
best_ptr->state  =  VISITED; 
delete_heap(  f rontier_heap,  best_ptr  ); 

return (  best_ptr  ); 
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tsearchx 

/* 

File:  tsearch.c 

Programmer:  g.b.  parker 

Environment :  any 

Language :  C 

Date:         9  July  92 

Revised: 

Comments:  This  is  a  multi-heuristic  search  that  takes  in  the  bias 

factors  in  the  form  of  an  eight  digit  hexadecimal  number. 
*/ 

♦include  "ga_search.h" 
♦include  <sys/time.h> 

int  heap_size; 


/*  ************** ***  search  *****************************************  */ 

float  search (  sx,  sy,  gx,  gy,  ind_chrom_f actor,  random_seed,  density,  node, 

ticks  ) 

int  sx,sy,gx,gy;     /*  position  in  density  array  for  start  &  goal  */ 

struct  f actor_struct  ind_chrom_f actor; 

int  random_seed; 

int  density[16]  [16]  ; 

struct  node_rec  *node [ 66] [66] ; 

float  *ticks; 

{ 

struct  node_rec  *current,  *previous,  *start,  *goal; 

struct  node_rec  *f rontier_heap [4096] ; 

long  sec,  usee; 
static  struct  timeval  *tvp; 
static  struct  timezone  *tzp; 
static  int  first  =  T; 

float  dist_traveled  =  0.0; 
int  k; 

if(  first  )  { 

tvp  =  (struct  timeval  *) malloc (sizeof (struct  timeval)); 

tzp  =  (struct  timezone  *)malloc (sizeof (struct  timezone)); 

first  =  F; 
} 
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srandom(random_seed) ;   /*   seed  the  random  generator  */ 

heap_size  =  0; 

make_array (density,  node); 

start  =  f ind_node (node,  sx,  sy) ; 
start->state  =  START; 
s t art ->dist_from_s tart  -  0.0; 
goal  =  f ind_node (node,  gx,  gy) ; 
goal->state  =  GOAL; 

for  (k=0;k<4096;k++) 

f rontier_heap[k]  =  NULL; 
current  =  start; 
previous  =  start; 

#ifdef  IRIS 

/*  initialize  the  IRIS  system  */ 

initialize (GATITLE) ; 
#endif 

gettimeofday (tvp,  tzp); 
sec  =  tvp->tv_sec; . 
usee  =  tvp->tv_usec; 

while  (  current  !=  goal)  { 

if (  adjacent (current ,  goal)  )  { 

dist_traveled  =  dist_traveled  +  update_dist (  current,  goal  ); 
current  =  goal; 
} 
else  { 

update_f rontier_list (  node,  current,  previous,  front ier_heap,  start, 

goal,  ind_chrom_factor  ) ; 
if (  heap_size  ==  0)  { 

print f("\nENDING  SEARCH  BEFORE  GOAL  -  no  more  frontiers- 
break; 
} 

previous  =  current; 
current  =  pick_best_f rontier (  node,  current,  front ier_heap,  goal, 

ind_chrom_f actor  ) ; 
dist_t raveled  =  dist_traveled  +  current->dist_f rom_current; 
} 
#ifdef  IRIS 

draw_terrain (node, start, goal, current, dist_traveled, ind_chrom_f actor) ; 
#endif 

)  /*  end  while  loop  */ 

gettimeofday (tvp,  tzp); 

*ticks  =  (float) (tvp->tv  sec  -  sec)  +  (tvp->tv_usec  -  usee) /1000000 . 0; 
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♦ifdef  IRIS 

sleep (  SLEEPTIME  ) ; 
#endif 

goal->state  =  GOAL; 
start->state  =  START; 

#ifdef  SUN 

/*  Print  to  standard  output  */ 

/*  not  normally  used,  but  optional  for  sun  */ 

/* 

printf ("\n") ; 

print_node (node) ; 

printf ("\n") ; 

printf ("\nDIST  =  %f",  dist_t raveled) ; 
*/ 
#endif 

return (  dist_traveled  ); 

} 
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tfrontier.c 

/* 

File:  tfrontier.c 

Programmer:  g.b.  parker 

Environment :  any 

Language :  C 

Date:  9  July  92 
Revised: 

Comments:  Maintenance  of  frontier  list 
*/ 

♦include  "ga_search.h" 
# include  "math.h" 


/*  *****************  update  frontier  list  ***************************  */ 
/*  Looks  two  away  from  the  current  node  to  update  stable  search 
characteristics  */ 

update_f rontier_list (  node,  c,  p,  f rontier_heap,  start,  goal,  factor  ) 

struct  node_rec  *node [66] [66] ; 

struct  node_rec  *c;         /*  current  */ 

struct  node_rec  *p;         /*  previous  */ 

struct  node_rec  *f rontier_heap [ 4096] ; 

struct  node_rec  *start; 

struct  node_rec  *goal; 

struct  f actor_struct  factor; 

{ 

int  xi,  yi,  base_xi,  base_yi,  top_xi,  top_yi; 

float  old_subtotal; 

base_xi  =  c->xi  ==  1  ?  1  :  c->xi  -  2 

base_yi  =  c->yi  ==  1  ?  1  :  c->yi  -  2 

top_xi  =  c->xi  ==  64  ?  64  :  c->xi  +  2 

top_yi  =  c->yi  ==  64  ?  64  :  c->yi  +  2 

update  adjacent  obstacles (  node,  c  ) ; 


94 


for  (xi=base_xi;xi<=top_xi;xi++)  { 

for  (yi=base_yi;yi<=top_yi;yi++)  { 

if  (  ( (node[xi] [yi]->state  ==  UNTOUCHED)  || 

(node [xi] [yi] ->state  ==  GOAL)  )   &&  adjacent (node [xi] [yi] , c)  ){ 
nodefxi] [yi]->state  =  FRONTIER; 
node [xi] [yi] ->predecessor  =  c; 

node [xi] [yi] ->dist_f rom_goal   =  update_dist (  node [xi] [yi] ,  goal  ); 
node [xi] [yi] ->dist_f rom_start   =  update_dist (  node [xi] [yi] ,  start  ); 
node [xi] [yi] ->subtotal  = 

compute_subtotal (  node [xi] [yi] ,  factor, 

calc_momentum(node [xi] [yi],c,p), 
update_crowd_sides (node, node [xi] [yi] ) , 

update_crowd_diag (node,  node [xi]  [yi] )  ); 
insert_heap(  f rontier_heap,  node[xi][yi]  ); 
) 

else  if  (  node [xi] [yi] ->state  ==  FRONTIER  )  { 
old_subtotal  =  node [xi] [yi] ->subtotal; 
node [xi] [yi] ->subtotal  = 

compute_subtotal (  node [xi] [yi] ,  factor, 

calc_momen turn  (node  [xi]  [yi]  ,  c,p)  , 
update_crowd_sides (node, node [xi] [yi] ) , 
update_crowd_diag (node, node [xi] [yi] )  ) / 
if  (  lequalf (old_subtotal, node [xi] [yi] ->subtotal)  )  { 

move_heap (  f rontier_heap,  node [xi] [yi] ->f rontier_index  )  ; 
} 
} 
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/*  *****************  update_adjacent_obstacles  **********************  */ 
/*  Records  for  future  use  which  adjacent  nodes  are  obstacles  */ 

update_adjacent_obstacles (  node,  c  ) 

struct  node_rec  *node [66] [66] ; 

struct  node_rec  *c;         /*  current  */ 

{ 

int  xi,  yi,  base_xi,  base_yi,  top_xi,  top_yi; 

base_xi  =  c->xi  -  1; 

base_yi  =  c->yi  -  1; 

top_xi  =  c->xi  +  1; 

top_yi  =  c->yi  +  1; 

for  (xi=base_xi;xi<=top_xi;xi++) 
for  (yi=base_yi ; yi<=top_yi ; yi++) 

if  (  nodefxi] [yi]->state  ==  OBSTACLE  ) 

node [xi] [yi] ->back_track_state  =  OBSTACLE; 

} 


/*  a****************  pick  best  frontier  *****************************  */ 
/*  finds  best  frontier  node,  returns  it  */ 

struct  node_rec  *pick_best_f rontier (node, current, f rontier_heap, goal, factor) 

struct  node_rec  *node [66]  [66]  ; 

struct  node_rec  *current; 

struct  node_rec  *f rontier_heap [4096] ; 

struct  node_rec  *goal; 

struct  factor_struct  factor; 

{ 

struct  node_rec  *best_ptr,  *q,  *qend,  *qreset; 

float  node_cost  =  BIG_NUMBER; 

float  norm  =  1.0;  /*  (current->dist_f rom_goal  /  16.0);    normalize  factor  */ 

float  lower_bound  =  f rontier_heap  [0] ->subtotal; 

float  upper_bound  =  BIG_NUMBER; 

int  xi,  yi,  k; 

float  steps; 

best_ptr  =  current; 
q  =  current; 
qend  =  current; 
qreset  =  current; 

current->back_track_state  =  VISITED; 
current->dist  from  current  =  0.0; 
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while  (  ( lower_bound  <  upper_bound)  &&  (q  !=  NULL)  )  { 

steps  =  q->dist_f rom_current  +  1.0; 
for(k=0;k<8;k++)  { 

if  (k==4) 

steps  =   q->dist_f rom_current  +  SQRT2; 

xi  =  gen_xi (  k,  q->xi  ) ; 

yi  =  gen_yi (  k,  q->yi  ); 

if  (  (  (nodefxi] [yi]->state  ==  VISITED)  || 
(nodefxi] [yi]->state  ==  FRONTIER)  || 
(nodefxi] [yi]->state  ==  START)  )  && 
nodefxi] [yi] ->back_track_state  ==  UNTOUCHED  )  { 
node [xi] [yi] ->dist_f rom_current  =  steps; 
node [xi] [yi] ->back_track_state  =  VISITED; 
node [xi] [yi] ->qreset  »  qreset; 
qreset  =  node [xi] [yi] ; 

if  (  (node[xi]  [yi]->state  ==  VISITED)  I  I 
(nodefxi] [yi]->state  ==  START)  )  { 
qend->qnext  =  node[xi] [yi] ; 
qend  =  node[xi] [yi] ; 
} 

else  {     /*  node[xi] [yi]->state  ==  FRONTIER  */ 
node_cost  =  node [xi] [yi] ->subtotal  + 

node [xi] [yi] ->dist_f rom_current  *  factor .current_dist  + 
calc_move_away (node [xi] [yi] , current, goal)  *  f actor .move_away; 
if  (node_cost  <  upper_bound)  { 
upper_bound  =  node_cost; 
best_ptr  =  node [xi] [yi] ; 
}  /*  end  if  */ 
}  /*  end  else  */ 
}  /*  end  if  */ 
}  /*  end  for  loop  */ 
q  =  q->qnext; 

lower_bound  =  f rontier_heap [0] ->subtotal  +  steps  *  factor .current_dist; 
}  /*  end  while  loop  */ 

reset_back_track_state (  qreset  ) ; 
best_ptr->state  =  VISITED; 
delete_heap(  f rontier_heap,  best_ptr  ); 

return (  best_ptr  ); 
} 
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/*  *****************  reset  back  track  state  *************************  */ 
/*  Resets  node  record  fields  used  to  perform  the  backtrack  search  */ 

reset_back_track_state (  qreset  ) 

struct  node_rec  *qreset; 

{ 

struct  node_rec  *temp; 

while  (qreset  !=  NULL)  { 
temp  =  qreset ->qreset; 
qreset->back_track_state  =  UNTOUCHED; 
qreset->qreset  =  NULL; 
qreset ->qnext  =  NULL; 
qreset  =  temp; 
1 
I 


/*  *****************  update  list  ************************************  */ 
/*  Euclidean  distance  between  input  nodes  */ 

float  update_dist (  nl,  n2  ) 
struct  node_rec  *nl; 
struct  node_rec  *n2; 
{ 

float  x  =  nl->x  -  n2->x; 

float  y  =  nl->y  -  n2->y; 

return (  sqrt (  x*x  +  y*y  )  ) ; 
} 


I  *  *****************  diff  int  ***************************************  */ 
/*  Absolute  difference  between  two  integers  */ 

dif f_int (  a,  b) 

int  a,  b; 

{ 

int  c  =  a  -  b; 

if  (  c  <  0  ) 

return (  -c  ) ; 
else 

return  (  c  ) ; 
} 
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'*  ***************** 


adjacent  ***************************************  */ 


I*   Returns  T  if  the  input  nodes  are  adjacent  */ 

adjacent (nl,n2) 
struct  node_rec  *nl; 
struct  node_rec  *n2; 

{ 

return  (  (dif f_int (nl->xi,n2->xi)  <  2)  &&  (dif f_int (nl->yi,n2->yi)  <  2)  ); 

} 


/*  *****************  update  crowd  sides  *****************************  */ 
/*  Counts  the  known  adjacent  horizontal/vertical  obstacles  to  the  frontier 
node  */ 


update_crowd_sides (node,  f) 
struct  node_rec  *node [ 66]  [ 66]  ; 
struct  node_rec  *f;   /*  frontier  */ 
{ 

int  s_count  =  0 ; 

/*  N  */ 

if  (  node [f->xi] [f->yi  +  1] ->back_track_state  ==  OBSTACLE  ) 

s_count  =  s_count  +  1; 
/*  E  */ 
if  (  node[f->xi  +  1] [f->yi] ->back_track_state  ==  OBSTACLE  ) 

s_count  =  s_count  +  1; 
/*  S  */ 
if(  node [f->xi] [f->yi  -  1] ->back_track_state  ==  OBSTACLE  ) 

s_count  =  s_count  +  1; 
/*  W  */ 
if(  node[f->xi  -  1] [f->yi] ->back_track_state  ==  OBSTACLE  ) 

s_count  =  s_count  +  1; 
return (  s_count  ) ; 
J 
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/*  *****************  update  crowd  diag  ******************************  */ 
/*  Counts  the  known  adjacent  diagonal  obstacles  to  the  frontier  node  */ 

upda t e_c ro wd_diag ( node ,  f ) 
struct  node_rec  *node [ 66] [ 66] ; 
struct  node_rec  *f;   /*  frontier  */ 
i 

int  d_count  =  0; 

/*  NE  */ 

if  (  node[f->xi  +  1] [f->yi  +  1] ->back_track_state  ==  OBSTACLE  ) 

d_count  =  d_count  +  1; 
/*  SE  */ 
if  (  node[f->xi  +  1] [f->yi  -  1] ->back_track_state  ==  OBSTACLE  ) 

d_count  =  d_count  +  1 ; 
/*  SW  */ 
if(  node[f->xi  -  1] [f->yi  -  1] ->back_track_state  ==  OBSTACLE  ) 

d_count  =  d_count  +  1; 
/*  nw  */ 

if(  node[f->xi  -  1] [f->yi  +  1] ->back_track_state  ==  OBSTACLE  ) 

d_count  =  d_count  +  1; 
return  (  d_count  ) ; 

) 


/*  *****************  calc  move  away  *********************************  */ 
/*  Determines  if  a  move  to  the  frontier  would  be  moving  away  from  the 
goal.   Each  axis  move  away  counts  as  two.  */ 

calc_move_away {  f,  c,  g) 
struct  node_rec  *f; 
struct  node_rec  *c; 
struct  node_rec  *g; 

{ 

int  ma_count  =  0 ; 

if(  diff_int (f->xi,g->xi)  >  dif f_int (c->xi, g->xi)  ) 

ma_count  =  2; 
if(  diff_int (f->yi,g->yi)  >  dif f_int (c->yi,g->yi)  ) 

ma_count  =  ma_count  +  2; 

return (  ma_count  ) / 
} 
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/*  *****************  calc  momentum  **********************************  */ 
/*  Returns  0  if  no  change  in  direction,  1  if  45  degree  change,  two  if  90 
degree  change,  and  three  if  135  degree  change  or  node  not  adjacent  */ 

calc_momentum(  f,  c,  p  ) 
struct  node_rec  *f; 
struct  node_rec  *c; 
struct  node_rec  *p; 
{ 

if(  adjacent (f, c)  &&  adjacent (c, p)  ) 

return  (  dif f_int (p->xi  -  c->xi,c->xi  -  f->xi)  +  dif f_int (p->yi  -  c->yi,c- 
>yi  -  f->yi)  ) ; 

else 

return (  3  ) ; 
} 


/*  *****************  compute  subtotal  *******************************  */ 
/*  Computes  the  frontier  nodes  subtotal  value  dependent  on  the  stable  heuris- 
tics */ 

float  compute_subtotal (  n,  factor,  m,  cs,  cd  ) 

struct  node_rec  *n;        /*  the  node  */ 

struct  factor_struct  factor; 

int  m;  /*  momentum  */ 

int  cs;  /*  crowding_sides  */ 

int  cd;  /*  crowding_diagonals  */ 

( 

return  (  n->dist_f rom_start  *  factor . start_dist  + 

n->dist_f rom_goal  *  factor .goal_dist  + 

cs  *  factor . crowd_sides  + 

cd  *  factor .crowd_diag  + 

m  *  factor .momentum  ); 
) 
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theap.c 

/* 

File:  theap.c 

Programmer:  g.b.  parker 

Environment:  any 

Language:  C 
Date:  9  July  92 

Revised: 

Comments:  Frontier  heap  functions 
*/ 

♦include  "ga_search .h" 

int  heap  size; 


/*  *****************  insert  heap  ************************************  */ 
/*  Inserts  a  node  into  the  frontier  heap  */ 

insert_heap(  fh,  n  ) 

struct  node_rec  *fh[4096];   /*  f rontier_heap  */ 

struct  node_rec  *n;         /*  node  to  insert  */ 

I 

n->state  =  FRONTIER; 

n->f rontier_index  =  heap_size; 

fh [heap_size]  =  n; 

heap_size  =  heap_size  +  1; 

move_heap (  fh,  heap_size-l  ); 
} 


/*  *****************  delete  heap  ************************************  */ 
/*  Deletes  a  node  from  the  frontier  heap  */ 

delete_heap(  fh,  n  ) 

struct  node_rec  *fh[4096];   /*  f rontier_heap  */ 

struct  node_rec  *n;         /*  node  to  delete  */ 

{ 

heap_size  =  heap_size  -  1; 

fh [n->f rontier_index]  =  fh [heap_size] ; 

fh [n->f rontier_index] ->f rontier_index  =  n->f rontier_index; 

n->state  =  VISITED; 

fh [heap_size]  =  NULL; 

if  (  n->f rontier_index  !=  heap_size  ) 
move_heap (  fh,  n->f rontier_index  ); 

n->f rontier_index  =  NA; 
} 


102 


/*  *****************  move  heap  **************************************  */ 
/*  Moves  a  node  in  the  frontier  heap  if  required;  dependent  on  the  value 
of  the  nodes  subtotal  field.  */ 

move_heap (  fh,  i  ) 

struct  node_rec  *fh[4096];   /*  f rontier_heap  */ 

int  i;  /*  index  of  node  to  possibly  move  */ 

{ 

int  parent  =  (i  -  1)  /  2; 

int  child  =  (  (2*i  +  l)  >=  heap_size)  ?  i  :  2*i  +  l; 

int  child2  =  (  (2*i+2)  >=  heap_size)  ?  i  :  2*i+2; 

if  (  (child2  !=  i)  &&  (fh [child2] ->subtotal  <  fh [child] ->subtotal)  ) 

child  =  child2; 
while  (fh[i] ->subtotal  <  fh [parent ] ->subtotal)  { 

swap(  fh,  i,  parent  ); 

i  =  parent; 

parent  =  (i  -  1)  /  2; 
} 
while  (fh[i]->subtotal  >  fh [child] ->subtotal)  { 

swap (  fh,  i,  child  ); 

i  =  child; 

child  =  ((2*i+l)  >=  heap_size)  ?  i  :  2*i+l; 

child2  -  ((2*i+2)  >=  heap_size)  ?  i  :  2*i+2; 

if  (  (child2  !=  i)  &&  (fh [child2] ->subtotal  <  fh [child] ->subtotal)  ) 
child  =  child2; 
} 
I 


/*  *****************  swap  *******************************************  *. 
/*  Swaps  nodes  in  the  frontier  heap  */ 

swap(  fh,  il,  i2  ) 

struct  node_rec  *fh[4096];   /*  f rontier_heap  */ 

int  il,  i2;  /*  indexes  of  nodes  to  swap  */ 

I 

struct  node_rec  *temp_ptr  =  fh[il]; 

fh [il] ->f rontier_index  =  i2; 
fh [i2] ->f rontier_index  =  il; 
fh[il]  =  fh[i2] ; 
fh[i2]  =  temp_ptr; 
I 
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evolve.c 

/* 


File:  evolve.c 

Programmer:  g.b.  parker 

Environment :  any 

Language :  C 

Date:  9  July  92 
Revised: 

Comments:  Performs  selection,  crossover,  and  mutation  on  input  population 
*/ 

#include  "ga_search.h" 


/*  *****************  evolve  *****************************************  */ 
/*  takes  in  a  population,  with  fitness  information  which  is  used  to  produce 
the  next  population  */ 

evolve (  individual,  rs  ) 

struct  individual_struct  *individual [32] ; 

int  rs;       /*  random_n umber  */ 

{ 

static  struct  individual_struct  *temp_ind [32] ; 

static  int  first  =  T;   /*  T  or  F  */ 

int  k; 

int  top  =  0; 

int  next_spot  =  2 ; 

int  mate_heap [31] ; 

int  even  =  T; 

if  (  first  )  { 

create_population (  temp_ind  ); 

first  =  F; 
} 

temp_ind[0] ->f it_sum  =  0.0; 
temp_ind[l] ->f it_sum  =  0.0; 


create_mate_heap (  mate_heap,  individual,  rs  ); 
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for  (k=0;k<32;k++)  { 

top  =  top  +  individual [k] ->fitness; 

while  (  (mate_heap[l]  <=  top)  &&  (mate_heap [0]  >  0)  )  ( 

temp_ind[next_3pot ] ->chrom. alleles  =  individual [k] ->chrom.  alleles; 
temp_ind [next_spot ] ->previous_index  =  k; 
if  (  even  )  { 

next_spot  =  next_spot  +2; 
even  =  next_spot  <  30  ?  T  :  F; 
} 
else 

next_spot  =  get_odd ( ) ; 
pop_mate_heap (mate_heap) ; 
} 
if  (  individual [k] ->fit_sum  >  temp_ind[0] ->f it_sum  )  { 

set_equal (  temp_ind[0] ,  temp_ind[l],  temp_ind [0] ->previous_index  ) ; 
set_equal (  individual [k] ,  temp_ind[0],  k  ); 
) 
else  if  (  (individual [k] ->fit_sum  >  temp_ind[l] ->f it_sum)  && 

(individual [k] ->chrom. alleles  !=  temp_ind[0] ->chrom. alleles)  ) 
set_equal(  individual [k] ,  temp_ind[l],  k  ) ; 
} 

crossover  (  individual,  temp_ind,  rs  )/ 
/*  mutate  done  in  crossover  */ 


/*  *****************  create  mate  heap  *******************************  */ 
/*  Creates  a  heap  of  integers  which  will  be  used  to  stochastically  choose 
individuals  for  reproduction  */ 

create_mate_heap (  mh,  ind,  rs  ) 

int  mh [31] ; 

struct  individual_struct  *ind[32]; 

int  rs;       /*  random_seed  */ 

I 

int  k; 

int  total_fit  =  0; 

s random (rs) ; 

mh[0]  =  0; 

for  (k=0;k<32;k++) 

total_fit  =  total_fit  +  ind[k] ->f itness; 

for  (k=l;k<31;k++) 

insert_mate_heap(  mh,  (random!)  %  total_fit)  ); 
I 
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/************************  crossover  functions  **************************/ 


/*  *****************  allele  crossover  *******************************  */ 
/*  Performs  crossover  of  the  chromosome  at  a  random  allele  position  */ 

allele_crossover (  ind,  temp_ind,  k,  cross_allele  ) 

struct  individual_struct  *ind[32]; 

struct  individual_struct  *temp_ind[32] ; 

int  k,  cross_allele; 

{ 

ind [k] ->chrom. alleles  =  temp_ind[k] ->chrom. alleles; 

ind[k+l] ->chrom. alleles  =  temp_ind [k+1] ->chrom. alleles; 

switch  (  cross_allele  ) 
{ 
case  0 : 

ind[k] ->chrom. f actor . start_dist  =  temp_ind[k+l] ->chrom. f actor . start_dist. 

ind [k+1] ->chrom. factor . start_dist  =  temp_ind[k] ->chrom. factor . start_dist, 
case  1: 

ind[k] ->chrom. factor .goal_dist  =  temp_ind[k+l] ->chrom. factor . goal_dist ; 

ind [k+1] ->chrom. factor .goal_dist  =  temp_ind [k] ->chrom. factor . goal_dist ; 
case  2 : 

ind[k] ->chrom. factor .cur rent_dist  = 

temp_ind[k+l] ->chrom. factor .cur rent_dist; 

ind [k+1] ->chrom. factor .current_dist  = 

temp_ind[k] ->chrom. factor .current_dist; 
case  3: 

ind[k] ->chrom. factor .crowd_sides  = 

temp_ind[k+l] ->chrom. f actor .crowd_sides; 

ind [k+1] ->chrom. factor . crowd_sides  = 

temp_ind [k] ->chrom. factor . crowd_sides; 
case  4: 

ind[k] ->chrom. factor .crowd_diag  =  temp_ind[k+l] ->chrom. factor .crowd_diaga 

ind [k+1] ->chrom. factor .crowd_diag  =  temp_ind[k] ->chrom. factor . crowd_diagj 
case  5: 

ind[k] ->chrom. f actor .move_a way  =  temp_ind [k+1] ->chrom. f actor .move_a way; 

ind [k+1] ->chrom. f actor .move_a way  =  temp_ind[k] ->chrom. f actor .move_a way; 
case  6: 

ind [k] ->chrom. factor .momentum  =  temp_ind [k+1] ->chrom. factor .momentum; 

ind [k+1] ->chrom. factor .momentum  =  temp_ind[k] ->chrom. factor .momentum; 
1 
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/*  *****************  bit  crossover  **********************************  */ 
/*  Performs  crossover  of  an  allele  at  a  random  bit  position  */ 

bit_crossover (  ind,  temp_ind,  k,  cross_allele  ) 
struct  individual_struct  *ind[32]; 
struct  individual_struct  *temp_ind [32] ; 
int  k,  cross_allele; 

{ 

int  cross_bit,  inv_cross_bit ; 

cross_bit  =  get_mask (  rand5 ( )  ) ; 
inv  cross  bit  =  cross  bit  *  MASK4; 


switch  (  cross_allele  ) 
{ 
case  1 : 

ind[k] ->chrom. factor . start_dist  = 

(  temp_ind[k] ->chrom. factor . start_dist  &  cross_bit  )    I 
(  temp_ind [k+1] ->chrom. factor . start_dist  &  inv_cross_bit  ); 
ind [k+1] ->chrom. f actor . start_dist  = 

(  temp_ind [k+1] ->chrom. factor . start_dist  &  cross_bit  )    I 
(  temp_ind [k] ->chrom. factor . start_dist  &  inv_cross_bit  )/ 
break; 
case  2 : 

ind[k] ->chrom. factor .goal_dist  = 

(  temp_ind [k] ->chrom. factor .goal_dist  &  cross_bit  )    I 
(  temp_ind [k+1] ->chrom. factor .goal_dist  &  inv_cross_bit  ); 
ind [k+1] ->chrom. factor .goal_dist  = 

(  temp_ind [k+1] ->chrom. factor .goal_dist  &  cross_bit  )    I 
(  temp_ind [k] ->chrom. factor .goal_dist  &  inv_cross_bit  ); 
break; 
case  3: 

ind [k] ->chrom. factor .current_dist  = 

(  temp_ind [k] ->chrom. factor .current_dist  &  cross_bit  )    I 
(  temp_ind [k+1] ->chrom. factor .current_dist  &  inv_cross_bit  ); 
ind [k+1] ->chrom. factor .current_dist  = 

(  temp_ind [k+1] ->chrom. factor .current_dist  &  cross_bit  )    I 
(  temp_ind[k] ->chrom. factor .current_dist  &  inv_cross_bit  ); 
break; 
case  4 : 

ind[k] ->chrom. factor . crowd_sides  = 

(  temp_ind[k] ->chrom. factor .crowd_sides  &  cross_bit  )     I 
(  temp_ind [k+1 ] ->chrom. factor .crowd_sides  &  inv_cross_bit  ); 
ind [k+1] ->chrom. factor . crowd_sides  = 

(  temp_ind [k+1] ->chrom. factor .crowd_sides  &  cross_bit  )     I 
(  temp_ind[k] ->chrom. factor .crowd_sides  &  inv_cross_bit  ); 
break; 
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case  5: 

ind[k] ->chrom. factor .crowd_diag  = 

(  temp_ind[k] ->chrom. factor. crowd_diag  &  cross_bit  )     I 
(  temp_ind[k+l] ->chrom. factor .crowd_diag  &  inv_cross_bit  ); 

ind [k+1] ->chrom. factor . crowd_diag  = 

(  temp_ind[k+l] ->chrom. factor .crowd_diag  &  cross_bit  )     I 
(  temp_ind[k] ->chrom. factor .crowd_diag  &  inv_cross_bit  ); 

break; 
case  6: 

ind [k] ->chrom. factor .move_away  = 

(  temp_ind [k] ->chrom. factor .move_away  &  cross_bit  )     I 
(  temp_ind [ k+1 ] ->chrom. factor .move_away  &  inv_cross_bit  ); 

ind[k+l] ->chrom. factor .move_away  = 

(  temp_ind[k+l] ->chrom. factor .move_away  &  cross_bit  )     I 
(  temp_ind[k] ->chrom. factor .move_away  &  inv_cross_bit  ); 

break; 
case  7 : 

ind [k] ->chrom. factor .momentum  = 

(  temp_ind [k] ->chrom. factor .momentum  &  cross_bit  )     I 
(  temp_ind [k+1] ->chrom. factor .momentum  &  inv_cross_bit  ); 

ind [k+1] ->chrom. factor .momentum  = 

(  temp_ind [k+1] ->chrom. factor .momentum  &  cross_bit  )     I 
(  temp_ind[k] ->chrom. factor .momentum  &  inv_cross_bit  ); 
1 


) 


I  *  *****************  crossover  **************************************  */ 
/*  The  main  function  */ 

crossover (  individual,  temp_ind,  rs  ) 
struct  individual_struct  * individual [32] ; 
struct  individual_struct  *temp_ind [32] ; 
int  rs;       /*  random_seed  */ 
{ 

int  k; 

int  cross_allele; 

s random (rs) ; 
cross_allele  =  rand8(); 

set_equal (  temp_ind [ 0 ] ,  individual [ 0 ] ,  temp_ind [ 0 ] ->previous_index  ) ; 
set_equal(  temp_ind[l],  individual [1] ,  temp_ind[l] ->previous_index  ); 

for(k=2;k<32;k=k+2)  { 

allele_crossover (  individual,  temp_ind,  k,  cross_allele  ); 
bit_crossover (  individual,  temp_ind,  k,  cross_allele  ); 
mutate  (  individual,  k,  rs  ) ; 

} 
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/*  *****************  get  mask  ***************************************  */ 
/*  called  by  bit_crossover  */ 

int  get_mask (  rn  ) 

int  rn;   /*  random  number  */ 

{ 

switch  (  rn  ) 

{ 

case  0 : 

return  MASKO; 
case  1 : 

return  MASK1; 
case  2: 

return  MASK2; 
case  3 : 

return  MASK3; 
case  4: 

return  MASK4; 
} 


/*  *****************  mutate  *****************************************  */ 
/*  Runs  through  each  bit  of  the  chromosome  determining  if  if  will  invert  */ 

mutate  (  ind,  k,  rs  ) 

struct  individual_struct  *ind[32]; 

int  k,  rs; 

( 

unsigned  int  mut_factorl  =  OxfffffffO; 

unsigned  int  mut_factor2  =  OxfffffffO; 

int  g; 

for  (g=0;g<28;g++)  { 

mut_factorl  =  (  mut_factorl  <<  1  )  +  one_if_mutate 0 ; 

mut_factor2  =  (  mut_factor2  <<  1  )  +  one_if_mutate () ; 
} 

ind [k] ->chrom. alleles  =  ind [k] ->chrom. alleles  A  mut_factorl; 
ind[k+l] ->chrom. alleles  =  ind [k+1] ->chrom. alleles  A  mut  factor2; 
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/*  *****************  one  if  mutate  **********************************  */ 
/*  Returns  1  if  mutation  is  to  take  place  at  the  present  bit  */ 

one_if_mutate ( ) 
{ 

if    (    randlOOOOO    <    PROB_BIT_MUTATE    ) 

return  (  1  )  ; 
else 

return (  0  ) ; 
} 


/*  *****************  get  odd  ****************************************  */ 
/*  Determines  placement  of  selected  individual  for  reproduction. 
Distributes  individuals  to  avoid  mating  of  like  chromosomes.  */ 

/*  Definitions  only  pertinent  to  this  function  */ 
#define  LOW  0 
#define  MED   1 
#define  HIGH  2 

get_odd ( ) 
{ 

static  int  next  =  LOW; 

static  int  base  =  1; 

switch  (  next  ) 
{ 
case  LOW: 

next  =  MED; 
base  =  base  +  2; 
return (  base  ) ; 
case  MED: 

next  =  HIGH; 
return  (  base  +  10  ); 
case  HIGH: 
next  =  LOW; 
if  (  base  ==  11  )  { 
base  =  1; 
return (  31  ) ; 
} 
else 

return (  base  +20  ); 
I 
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/*  *****************  set-equal  **************************************  */ 
/*  Sets  one  individual  equal  to  another  */ 

set_equal (  f rom_ind,  to_ind,  k  ) 

struct  individual_struct  *from_ind,  *to_ind; 

int  k; 

{ 

to_ind->chrom. alleles  =  f rom_ind->chrom. alleles; 

to_ind->f itness  =  f rom_ind->f itness; 

to_ind->f it_sum  =  f rom_ind->f it_sum; 

to_ind->previous_index  =  k; 
} 
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eheap.c 

/* 

File:  eheap.c 

Programmer:  g.b.  parker 

Environment:  any 

Language :  C 
Date:  9  July  92 

Revised: 

Comments:  Frontier  heap  functions 
*/ 

♦include  "ga_search.h" 


/*  *****************  insert  mate  heap  *******************************  */ 

insert_mate_heap (  mh,  num  ) 

int  mh[31];   /*  mate_heap  */ 

int  num;  /*  number  to  insert  */ 

{ 

mh[0]  =  mh[0]  +  1; 

mh [  mh[0]  ]  =  num; 

move_mate_heap (  mh,  mh[0]  ); 
} 


/*  *****************  poo  mate  heap  **********************************  */ 
/*  Removes  top  of  mate  heap  */ 

pop_mate_heap (  mh  ) 

int  mh[31];   /*  mate_heap  */ 

{ 

mh[l]  =  mh[  mh [ 0 ]  ]  ; 

mh[  mh[0]  ]  =  0; 

mh[0]  =  mh[0]  -  1; 

move_mate_heap (  mh,  1  ); 
} 
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/*  *****************  move  mate  heap  *********************************  */ 
/*  Readjusts  heap  after  addition /removal  of  one  of  its  members  */ 

move_mate_heap (  mh,  i  ) 

int  mh[31];     /*  mate_heap  */ 

int  i;  /*  index  of  num  to  possibly  move  */ 

t 

int  parent  =i==l?l:i/2; 

int  child  =  (  (2*i)  >  mh[0])  ?  i  :  2*i; 

int  child2  =  (  (2*i+l)  >  mh[0])  ?  i  :  2*i+l; 

if  (  (child2  !=  i)  &&  (mh[child2]  <  mh[child])  ) 
child  =  child2; 

while  (mh[i]  <  mh [parent])  { 

swap_num(  mh,  i,  parent  ) ; 

i  =  parent; 

parent  =i==l?l:i/2; 
} 

while  (mh[i]  >  mh [child])  { 
swap_num(  mh,  i,  child  ) ; 
i  =  child; 

child  =  ((2*i)  >  mh[0])  ?  i  :  2*i; 
child2  =  ((2*i+l)  >  mh[0])  ?  i  :  2*i+l; 
if  (  (child2  !=  i)  &&  (mh[child2]  <  mh[child])  ) 
child  =  child2; 
} 
} 


/*  *****************  swap  num  ***************************************  */ 
/*  Swaps  positions  of  two  members  of  the  mate_heap  */ 

swap_num(  mh,  il,  i2  ) 

int  mh[31];   /*  mate_heap  */ 

int  il,  i2;   /*  indexes  of  numbers  to  swap  */ 

{ 

int  temp_num; 

temp_num  =  mh[il]; 
mh[il]    =  mh[i2] ; 
mh[i2]    =   temp_num; 
} 


113 


tmiscc 

/* 

File:  tmiscc 

Programmer:  g.b.  parker 

Environment :  any 

Language :  C 

Date:  6  apr  92 
Revised: 

Comments:  Miscellaneous  functions 
*/ 

♦include  "ga_search .h" 


a****************************************  * , 


/*  *****************  qen  xi 

/*  Generates  an  integer  value  dependent  on  the  input  xi  and  k. 

Used  with  gen_yi  to  generate  all  adjacent  nodes  to  (xi,yi) .  */ 

int  gen_xi (  k,  xi  ) 
int  k; 
int  xi ; 
{ 

switch  (k) 
{ 

case  0  : 

return (  xi  ) ; 
case  1: 

return (  xi+1  ) ; 
case  2  : 

return (  xi  ) ; 
case  3: 

return (  xi-1  ) ; 
case  4: 

return (  xi+1  ) ; 
case  5: 

return (  xi+1  ) ; 
case  6: 

return (  xi-1  ) ; 
case  7 : 

return (  xi-1  ) ; 
} 
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/*  *****************  gen  yi  *****************************************  */ 
/*  Generates  an  integer  value  dependent  on  the  input  yi  and  k. 
Used  with  gen_xi  to  generate  all  adjacent  nodes  to  (xi,yi) .  */ 

int  gen_yi (  k,  yi  ) 
int  k; 
int  yi ; 
{ 

switch  (k) 

{ 

case  0 : 

return (  yi+1  ) / 
case  1: 

return (  yi  ) ; 
case  2 : 

return (  yi-1  ) ; 
case  3: 

return  (  yi  ) ; 
case  4 : 

return (  yi+1  ) ; 
case  5 : 

return  (  yi-1  ) ; 
case  6 : 

return  (  yi-1  ) ; 
case  7 : 

return (  yi+1  ) ; 
} 
} 

/*  *****************  ecrualf  *****************************************  */ 
/*  Checks  if  two  floats  are  equal  (within  0.0001)  */ 

int  equalf (  x,  y  ) 

float  x,  y; 

{ 

if  (  ( (x-y)  <  -0.0001)  ||  ( (x-y)  >  0.0001)  ) 

return (  F  ) ; 
else 

return  (  T  ) ; 
} 
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/*  *****************  show  least  nodes  *******************************  */ 
/*  Sets  the  state  field  to  x  for  all  nodes  in  the  shortest  path  */ 
/*  Not  currently  used,  but  available  for  graphics  */ 

show_least_nodes (  node,  g  ) 
struct  node_rec  *node [ 66] [ 66] ; 
struct  node_rec  *g;   /*  goal  */ 
{ 

struct  node_rec  *best_ptr; 

float  best; 

int  xi,yi,k; 

while  ( (g->dist_f rom_start  >  0.0)  &&  (g->dist_f rom_start  <  10000.0))  { 
best  =  BIG_NUMBER; 
for(k=0;k<8;k++)  { 

xi  =  gen_xi (  k,  g->xi  ) ; 
yi  =  gen_yi (  k,  g->yi  ) ; 
if(  (  node [xi] [yi] ->state  !=  OBSTACLE)  && 

(  node [xi] [yi] ->dist_f rom_start  <  best  )  )  { 
best  =  node [xi] [yi] ->dist_f rom_start ; 
best_ptr  =  node [xi] [yi] ; 
} 
} 

g  =  best_ptr; 
g->state  =  X; 
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tdisplay.c 

/* 

File:  tdisplay.c 

Programmer:  g.b.  parker 

Environment :  any 

Language:  C 

Date:  9  July  92 

Revised: 

Comments:  Functions  called  by  all  searches  to  display  the  search  on  the 
IRIS.   This  file  should  not  be  excluded  from  Makefile  if  compiled  on  the  SUN, 
*/ 

♦include  "ga_search.h" 
♦include  <gl.h> 
♦include  <device.h> 


/*  *****************  initialize  *************************************  */ 
/*  Initializes  graphics  systems  for  output  */ 

initialize (title) 
char  title [33] ; 
{ 

/*  set  up  a  preferred  size  and  location  for  the  window  */ 
pref size (XMAXSCREEN+1, YMAXSCREEN+1-25  6) ; 
prefposition(0, 980, 0, 980) ; 

/*  open  a  window  for  the  program  */ 
winopen ("search" ) ; 

/*  put  a  title  on  the  window  */ 
wintitle (title) ; 

/*  put  the  machine  into  double  buffer  mode  */ 
doublebuf f er  ( ) ; 

/*  set  RGB  mode  for  color  */ 
RGBmode ( ) ; 

/*  configure  the  IRIS  (means  use  the  above  command  settings)  */ 
gconf  ig()  ; 

/*  queue  the  redraw  device  */ 
qdevice (REDRAW) ; 

/*  queue  buttons  needed  */ 
qdevice (BUT6) ;   /*  ESC  */ 
qdevice (BUT50) ;   /*  enter  */ 
qdevice (BUT4) ;   /*  right  shift  */ 
qdevice (BUT7 3) ;  /*  down  arrow  */ 

/*  set  the  world  coordinate  system  */ 
ortho2  (-1. 0, 66. 0,-1.0, 66.0)  ; 
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/*  *****************  draw  terrain  ***********************************  */ 
/*  Called  by  searches  to  draw  the  node  array  */ 

draw_terrain (node, start, goal, current, dist, chrom) 

struct  node_rec  *node [66] [66] ; 

struct  node_rec  *start, *goal, *current ; 

float  dist;   /*  dist_traveled  */ 

struct  factor_struct  chrom;  /*  ind_chrom_f actor  */ 

( 

short  value; 

static  int  cont  =  T; 

int  mmouse  =  F; 

int  first  =  T; 

int  do_print; 

if(  adjacent (start, current)  II  adjacent (goal, current )  ) 
cont  =  T; 

while (  (mmouse  I  I  first)  &&  cont  )  { 
do_print  =  F; 
draw_grid () ; 
draw_nodes (node, start, goal, current) ; 

while  (  qtestO  ) 

switch (  qread(Svalue)  ) 
( 
case  BUT6:       /*  "ECS"  to  terminate  display  for  that  search  */ 

cont  =  F; 

break; 
case  BUT50:      /*  "return"  to  halt  display  */ 

mmouse  =  T; 

break; 
case  BUT4:       /*  "shift"  to  continue  display  */ 

mmouse  =  F; 

break; 
case  BUT73:      /*  "down  arrow"  to  print  node  info  to  standard  output 
*/ 

do_print  =  T;         /*  node  is  selected  by  mouse  position  */ 

break; 
default : 

break; 
} 
show_mouse (node, dist , chrom, do_print) ; 

swapbuf fers () ;      /*  change  the  buffers  ...  */ 
first  =  F; 

} 
} 
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/*  *****************  show  mouse  *************************************  */ 
/*  Shows  mouse  position  and  prints  node  info  if  selected  */ 

show_mouse (node, dist, c,do_print) 

struct  node_rec  *node [ 66] [66] ; 

float  dist; 

struct  f actor_struct  c;  /*  ind_chrom_f actor  */ 

int  do_print; 

( 

int  mx_pix  =  getvaluator (MOUSEX) ; 

int  my_pix  =  getvaluator (MOUSEY) ; 

int  mx  =  ((67  *  mx_pix) /980)  -  1; 
int  my  =  ((67  *  my_pix)/980)  -  1; 

if  (  mx  >  65  ) 

mx  =  65; 
if  (my  >  65  ) 

my  =  65; 

RGBcolor  (0, 0,0); 

square (  node [mx] [my] ->x,  node [mx] [my] ->y,  0.35  ); 

if (  do_print  )  ( 

printf ("\n  dist=%f",  dist); 

printf ("\n(%d, %d) \n%d  state\n%d  btstate\n%f  subtotal\n%f  %d  start 
\n%f  %d  goal\n%f  %d  current\n%d  frontier", 
nodefmx] [my]->xi,  node [mx] [my]->yi,  node [mx] [my] ->state, 
nodefmx] [my] ->back_track_state,  node [mx] [my] ->subtotal, 
nodefmx] [my] ->dist_f rom_start,  c . start_dist, 
nodefmx] [my] ->dist_f rom_goal,  c.goal_dist, 
node [mx] [my] ->dist_f rom_current,  c .current_dist, 
node [mx] [my] ->f rontier_index) ; 
printf ("\n%f  from  below",  node [mx] [my] ->subtotal  - 

(node [mx] [my] ->dist_f rom_start  *  c.start_dist  + 
node [mx] [my] ->dist_f rom_goal  *  c.goal_dist)  ); 
printf ("\ncs=%d,cd=%d,ma=%d,m=%d\n", 

c .crowd_sides, c . crowd_diag, c .move_away, c .momentum) ; 
} 
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/*  *****************  draw  nodes  *************************************  */ 
/*  Draws  node  info,  can  be  changed  for  color  */ 

draw_nodes (node, start,  goal,  current) 
struct  node_rec  *node [66]  [66]  ; 
struct  node_rec  *start, *goal, *current ; 
{ 

int  xi,  yi; 

for (xi=0;xi<=65;xi++) 
for(yi=0;yi<=65;yi++)  { 

switch (  node [xi] [yi] ->state  ) 
{ 
case  OBSTACLE: 

RGBcolor (0, 0,0) ; 

squaref(  node [xi] [yi] ->x,  node [xi] [yi] ->y,  0.5  ); 
break; 
case  VISITED: 

/*  RGBcolor (0,0,  255)  ;  */ 

circf(  node [xi] [yi] ->x,  node [xi] [yi] ->y,  0.1  ); 
circ(  node [xi] [yi] ->x,  node [xi] [yi] ->y,  0.3  ); 
break; 
case  FRONTIER: 

/*  RGBcolor(0, 255,0)  ;  */ 

circ(  node [xi] [yi] ->x,  node [xi] [yi] ->y,  0.3  ); 
break; 
} 
1 

/*  RGBcolor(255, 0,255)  ;  */ 

circf (  current->x,  current->y,  0.3  ); 

/*  RGBcolor (255,  0,  0)  ;  */ 

circf (  start->x,  start->y,  0.4  ); 

circf (  goal->x,  goal->y,  0.4  ); 

) 
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/*  *****************  draw  grid  **************************************  */ 
/*  Draws  the  cross  lines  for  the  grid  */ 

draw_grid ( ) 
{ 

int  i  ; 

float  fi; 

/*  draw  the  background  color  */ 
RGBcolor(255,255,255)  ; 
clear  ()  ; 

RGBcolor (0, 0,0); 

for  (i=0;i<=66;i++)  { 
fi  =  (float) (i-0.5) ; 
move2 (-0.5,fi) ; 
draw2 (65.5, fi) ; 
move2 (f i, -0  .  5)  ; 
draw2 (f i, 65.5) ; 
} 
} 

/*  *****************  scruaref  ****************************************  */ 
/*  display  filled  square  for  2D  displays  */ 


void  squaref (xc, yc,d) 

float  xc,yc;   /*  center  point  of  square  */ 

float  d;       /*  half  of  side  length  */ 

{ 

rectf (xc-d, yc-d, xc+d, yc+d) ; 
I 


/*  *****************  sq-uare  *****************************************  */ 

/*  display  square  for  2D  displays  */ 

void  square (xc,yc,d) 

float  xc,yc;   /*  center  point  of  square  */ 

float  d;      /*  half  of  sidelength  */ 

I 

rect (xc-d, yc-d, xc+d, yc+d) ; 
I 
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tprint.c 

/* 

File:  tprint.c 

Programmer:  g.b.  parker 

Environment:  any 

Language :  C 

Date:  9  July  92 

Revised: 

Comments:  Prints  to  standard  output 
*/ 

♦include  "ga_search.h" 

/*  *****************  print  density  **********************************  */ 
/*  Prints  density  terrain  to  standard  output  */ 

print_density (  density  ) 
int  density [16] [16] ; 
I 

int  i,  j; 

for  (  j  =  15;  j>=0;  j—  )  {• 
print f ("\n") ; 
for  (  i=0;  i<=15;  i++  ) 

printf("%x  ",  density [i] [j] ) ; 

} 
} 


/*  *****************  print  population  *******************************  */ 
/*  Prints  the  population  to  standard  output  */ 

print_population (  i  ) 

struct  individual_struct  *i[32]; 

{ 

int  k; 

for  (k=0;k<32;k++) 

printf("\n  %d  %x  %d  %f  %d",  k,  i [k] ->chrom. alleles,  i [k] ->f itness,  i [k] 
>fit_sum,  i [k] ->previous_index) / 

} 
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/*  *****************  print  node  *************************************  */ 
/*  Prints  node  terrain  to  standard  output  */ 

print_node (  node  ) 

struct  node_rec  *node [66] [66] ; 

{ 

int  i,  j; 

for  (  j  =  65;  j>=0;  j—  )  { 
printf  C\n") ; 
for  (  i=0;  i<=65;  i++  ) 

switch  (node [i] [ j] ->state) 
{ 
case  UNTOUCHED: 

printf (".") ; 

break; 
case  OBSTACLE: 

printf ("#") ; 

break; 
case  VISITED: 

printf ("o") ; 

break; 
case  FRONTIER: 

printf ("f")  ; 

break; 
case  START: 

printf ("S") ; 

break; 
case  GOAL: 

printf ("G") ; 

break; 
case  CURRENT: 

printf ("O") ; 

break; 
case  SHORTEST: 

printf ("s") ; 

break; 
case  X: 

printf ("x") ; 

break; 
} 
} 
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